Emergence of helicity ±2 modes (gravitons) from qubit models 
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Gauge symmetry is commonly regarded as one of tlie founding principles of nature. But recent 
studies of topological/quantum order suggest that gauge symmetry can emerge as a low energy 
property of a qubit model that have no gauge symmetry at all at lattice scale. This suggests that 
gauge symmetry may not be a founding principle, but merely a property of a quantum ground state 
with long range quantum entanglements. 

The general equivalence principle and the associated diffeomorphism gauge symmetry are also 
regarded as a founding principles of nature. So one may wonder, can diffeomorphism gauge sym- 
metry also emerge as a low energy property of certain topological/quantum order in a qubit model? 
In this paper, we designed qubit models (or quantum spin models) on 3D lattice and showed that, 
for the first time, gapless helicity ±2 excitations (ie the gravitons) can emerge as the only low 
energy excitation in our models. We showed that the emergence of gapless helicity ±2 excita- 
tions in our models leads to the emergence of, at least, linearized diffeomorphism gauge symmetry 
h^v — >■ hfiv + d^fv + dvfii at low energies. 

In the first qubit model (called the L-type model), we show that helicity ±2 gapless excitations 
appear as the only type of low energy excitations using a reliable semiclassical approach. The 
dispersion of the gapless helicity ±2 is found to be efc oc |fcj^. The appearance of the gapless helicity 
±2 modes suggests that the ground state of the qubit model is a new state of matter. In the 
second model (called the N-type model) the collective modes are strongly interacting and there is 
no reliable approach to understand its low energy dynamics. Using a spin-wave/quantum-freeze 
approach (which is shown to reproduce the correct emergent U{1) gauge theory in a quantum rotor 
model), we argue that the second model may contain helicity ±2 gapless excitations as the only 
type of low energy excitations with a linear dispersion cj oc fc. 

Both models have emergent low energy diffeomorphism gauge symmetry which leads to the as- 
sociated gapless helicity ±2 excitations. We believe that those properties are topologically robust: 
any translation invariant perturbations cannot break the emergent diffeomorphism gauge symmetry 
and cannot generate a gap for those helicity ±2 excitations. Our results shed light on the quest 
to find a quantum theory of gravity - a quantum model with a finite cutoff whose ground state 
supports gravitons (ie helicity ±2 gapless excitations with a linear ej, oc |fc| dispersion). 



Contents 

I. Emergence approach 

A. Discoveries and unifications 

B. Seven wonders of universe 

C. Locality principle and unification 

D. Two approaches 

E. Originate from organization 

II. The strategy of our approach 

A. The rule of game 

B. A brief outline of our approach 

C. Some previous approaches 

III. Review of U(l) gauge theory 

A. Removing helicity mode though Gauss 
constraint 

B. Try to gap helicity mode though energy 
penalty 

C. Another failed attempt 

D. Gapping helicity modes on lattice with 
discretized £i-j (or campactified Oij) 

E. Low energy collective modes of the rotor model 
through spin- wave/quantum- freeze approach 

IV. Generalization to symmetric tensor field 

A. A field theory of symmetric tensor 



B. diE^-' = constraint and helicity and ±2 
quadratic modes 10 

C. Additional R" = constraint to remove the 
second helicity mode 10 

D. Nonlocal gauge invariant Hamiltonian 11 



V. 



The lattice models and imposing constraints 
through energy penalties 

A. Putting the constraints on lattice 

B. Compactification 

C. Low energy Hilbert space 



VI. 



The L-type lattice model 

A. Putting the R^^ R'^ term on lattice 

B. Putting the CjCj term on lattice 

C. Low energy collective modes - 
spin- wave / quantum- freeze approach 

1. Classical spin waves 

2. The classical helicity ±2 modes 

3. The quantum freeze and the gapping of 
helicity and ±1 modes 



VII. The N-type lattice model 

A. Putting OijR^-' on lattice 

B. Putting £'^£'^ - on lattice 

C. The N-type lattice model 

D. Low energy collective modes - 
spin- wave / quantum- freeze approach 



12 
12 
13 
14 

15 
16 
16 

17 
17 
17 

17 

18 
18 
18 
18 

19 



2 



1. Classical spin waves 19 

2. The classical helicity ±2 modes 20 

3. The quantum freeze and the gapping of 
helicity and ±1 modes 20 

E. Difficulties with the N-type lattice model 20 

VIII. Topological defects 21 

IX. Conclusion and discussion 21 

A. Low energy collective modes for L-type 

lattice model — equation-of-motion approach 23 

References 25 



I. EMERGENCE APPROACH 
A. Discoveries and unifications 



• The development of theoretical physics is driven by 
the desire to understand everything from a single 
or very few origins. 

Although the down pull by the earth was realized even 
before human civilization, such a phenomena did not 
arose any curiosity. After Galileo and Kepler found that 
planets move in a certain particular way described by a 
mathematical formula, people started to wonder: why 
planets move in such a peculiar and precise way. This 
motivated Newton to develop his theory of gravity. New- 
ton's theory not only explains the planets motion, it also 
explains the down-pull that we feel on earth. It uni- 
fies the two seemingly unrelated phenomena. Later we 
discovered that two other seemingly unrelated phenom- 
ena, electricity and magnetism, can generate each other. 
Our curiosity about the electricity and magnetism leads 
to another giant leap in science, which is summarized 
by Maxwell equations. Maxwell theory unifies electricity 
and magnetism and reveals that light is merely an elec- 
tromagnetic wave. We gain a much deeper understanding 
of light, which is so familiar and yet so unexpectedly rich 
and complex in its internal structure. Newton's theory of 
gravity and Maxwell's theory of light illustrate how sci- 
ence develops: it develops through the cycles of discovery, 
unification, more discovery, more unification. 



B. Seven wonders of universe 



• The current physical understanding of our world is 
build on seven basic assumptions. 

• Field theory and geometry played the central role in 
formulating the fundamental theory of nature. 

• This paper is a non-geometric approach to gravity 
(and other wonders). 



Modern science has made many more discoveries and 
has also unified many seemingly unrelated discoveries 
into a few simple structures. Those simple structures 
are so beautiful and we regard them as wonders of our 
universe. They are also very mysteries since we do not 
understand where do they come from and why do they 
have to be the way they are. At moment, the most fun- 
damental mysteries and/or wonders in our universe can 
be summarized by the following short list: 

(1) Locality. 

(2) Identical particles. 

(3) Gauge interactions.^"'^ 

(4) Fermi statistics. ^'^ 

(5) Chiral fermions.^'^" 

(6) Lorentz invariance."'^"'^ 

(7) Gravity 12 

In the current physical theory of nature, we take the 
above properties for granted and do not ask where do 
they come from. We put those wonderful properties into 
our theory by hand, for example, by introducing one field 
for each kind of interactions or elementary particles. 

Here, we would like to question where those wonder- 
ful and mysterious properties come from. Following the 
trend of science history, we wish to have a single unified 
understanding of all of the above mysteries. Or more pre- 
cisely, we wish that we can start from a single structure 
to obtain all of the above wonderful properties. 

Historically, the first attempt to unify different forces 
is motivated by Einstein general relativity which view 
gravity as an distortion of space. In 1918, Weyl pro- 
posed that the unit that we used to measure physical 
quantities is relative and is defined only locally. A dis- 
tortion of the unit system can be described by a vec- 
tor field which is called gauge field. Weyl proposed that 
such a vector field (the gauge field) is the vector potential 
that describes the electromagnetism. Although the above 
particular proposal turns out to be incorrect, the Weyl's 
idea is correct. In 1925, the complex quantum amplitude 
was discovered. If we assume the complex phase is rela- 
tive, then a distortion of unit system that measure local 
complex phase can also be described by a vector field. 
Such a vector field is indeed the vector potential that 
describes the electromagnetism. This leads to a unified 
way to understand gravity and electromagnetism: grav- 
ity arises from the relativity of spacial directions at dif- 
ferent spatial points, while electromagnetism arises from 
the relativity of complex quantum phases at different spa- 
tial points. Nordstrom, Moglichkeit, Kaluza, and Klein 
further showed that both gravity and electromagnetism 
can be understood as a distortion of space provided that 
we think the space as five dimensional with one dimen- 
sion compactified into a small circle. ^^^^^ Since that time, 
the geometric understanding of gravity and electromag- 
netism have dominated theoretical physics. 

In this paper, 1® we will take a position that the ge- 
ometric understanding is not good enough and will try 
to advocate a very different non-geometric understand- 
ing of gravity and electromagnetism. Why the geometric 
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understanding is not good enough? First the geometric 
understanding is not self-consistent. The consideration 
based quantum mechanics and Einstein gravity indicates 
that two points separated by a distance less than the 
Planck length cannot exist as a physical reality. Thus 
the foundation of the geometric approach - manifold - 
simply docs not exist in our universe. This suggests that 
geometry is an emergent phenomenon that appears only 
at long distances. So we cannot use geometry and mani- 
fold as a foundation to understand fundamental physical 
problems. 

Maxwell theory of light and Einstein theory of gravity 
are built on top of geometry. They fail to answer what 
is the origin geometry. In other words, Maxwell theory 
and Einstein theory predict light waves and gravitational 
waves. But the theories fail to tell us what is waving? 
Also, Maxwell theory and Einstein theory fail to answer 
what is the origin of the other wonders of our universe, 
such as Fermi statistics. 



C. Locality principle and unification 

• Quantum theory and locality suggest a very natural 
fundamental build block - qubit. But, can qubits 
unify everything? 

In this paper we try to address the above question by 
trying to find a single microscopic origin for both light 
and gravity. We hope the microscopic origin will allow 
us to gain a deeper and a unified understanding of grav- 
ity and light. We also hope that a deeper and unified 
understanding of gravity and light will lead to a unified 
understanding of other deep mysteries of the universe, 
such as the origin of all elementary particles. 

At first sight, it appears that such a goal is too am- 
bitious to be practical. However, recent progresses sug- 
gest that such an ambitious goal may be achievable. It 
was shown that, starting from a single origin - a lo- 
cal bosonic model (which is also called a quantum spin 
model, or more precisely a qubit model), the first four 
of the seven wonders can emerge naturally at low ener- 
gies, if the ground state of the spin model is described 
by string-net condensation. Thus we may say that 
the string-net condensation in a qubit model provides a 
unified origin for the four wonders: identical particles, 
gauge interactions, Fermi statistics and near massless- 
ness of the fermions. Two seemingly unrelated proper- 
ties, Gauge interaction and Fermi statistics, arc unified 
under the string-net picture. It is amazing to see that us- 
ing emergence and the string-net picture to understand 
the origin of light leads to a understanding of the origin 
of Fermi statistics. 

Knowing the above result, the goal of this paper can 
be stated more precisely: in addition to the four emer- 
gent properties mentioned above, we want to understand 
merely one more emergent property - the emergence of 
gravity, from a local qubit model. In this paper, we 



want to construct a local qubit model where gravitational 
waves and gravitons emerge as the low energy excitations 
of the model. If successful, such a local qubit model will 
actually represent a quantum theory of gravity. So con- 
structing such a local qubit model solves a long stand- 
ing problem of putting quantum mechanics and gravity 
together. We hope in the future, we can construct a 
local qubit model where light, gravity, and nearly mass- 
less fermions emerge at low energies. Such a local qubit 
model will unify five of the seven deep mysteries in our 
list. 

The belief that all the wonderful phenomena of our uni- 
verse (such as gauge interaction, Fermi statistics, gravi- 
tons, or even superstring theory) emerge from a lowly 
qubit model is called locality principle.^ Using a local 
qubit model as a underlying structure to understand the 
deep mysteries of our universe^ ''■^^ represents a departure 
from the traditional approach to understand our world by 
dividing things into smaller parts (which will be called di- 
vision approach) . It also represents a departure from the 
geometric way to understand gravity and gauge interac- 
tions. 



D. Two approaches 

• In the traditional division approach, people are look- 
ing for the fundamental build block of matter, by 
dividing particles into smaller and smaller pieces. 

• In the emergence approach, we are looking for the 
fundamental build block of space itself and assume 
that such building block is qubit. The elementary 
particles are viewed as excitations of qubits that 
from the space. 

In the division approach, we try to understand various 
things by divide them into smaller and smaller parts. If 
we assume the division has to end at a certain level, then 
we conclude that all things are formed by the parts that 
cannot be divided further. The indivisible parts are call 
the elementary particles. So in the division approach, 
we view everything in our world as made of some sim- 
ple beautiful building blocks, the elementary particles. 
A deeper understanding is gained if we find some ele- 
mentary particles are not actually elementary and are 
formed by even smaller objects. A large part of science 
is devoted in finding those smaller and smaller objects, 
as represented by the discoveries of atoms, electrons and 
protons, and then quarks. 

However, the division approach that we followed in last 
150 years may not represent a right direction. For exam- 
ple, phonons in a solid is as particle-like as any other 
elementary particles at low energies. But if we look at 
phonons closely, we do not see smaller parts that form a 
phonon. We see the atoms that fill the entire space. The 
phonons are not formed by those atoms, the phonons are 
simply collective motions of those atoms. 
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This leads us to wonder that maybe photons, electrons, 
gravitons, etc , are also collective motions of a certain un- 
derlying structure that fill the entire space. They may 
not have smaller parts. Looking for the smaller parts of 
photons, electrons, and gravitons to gain a deeper under- 
standing of those elementary particles may not be a right 
approach. 

In this paper, we will use a different approach, emer- 
gence approach, to gain a deeper understanding of ele- 
mentary particles. In the emergence approach used in 
this paper, we view space as a collection of qubits. The 
empty space (the vacuum) corresponds to the ground 
state of the qubits, and the elementary particles (the 
matter formed by them) correspond to the excitations 
of the qubits. 

We know that elementary particles in our would can 
have quite different properties. Can excitations of qubits 
have those rich properties? Due to the particle-wave du- 
ality in quantum theory, particles and waves are the same 
thing. In the emergence approach we will try to under- 
stand the nature of waves in order to understand the 
nature of elementary particles. 

E. Originate from organization 

• In the emergence approach, the particular prop- 
erties of elementary particles do not come from 
the particular choice of building blocks (there is 
only one choice in the emergence approach ~ the 
qubit), but come from how qubits are organized in 
the ground state. 

To see how the emergence approach leads to a deeper 
understanding of elementary particles, let us start with 
Euler equation dfP — v^d^p = that describe waves in 
a liquid. We want to ask what is the microscopic origin 
of Eulcr equation? The answer is obtained only after the 
discovery of atoms. We find that the waves in a liquid 
are collective motions of the atoms and Euler equation 
describes the dynamics of the those collective motions. 
However, how atoms are organized in the ground state 
plays an important role here. Only when the atoms are 
organized into a liquid state (or more precisely, a boson 
condensed state for the bosonic atoms), are the collec- 
tive motions of those bosonic atoms described by Euler 
equation. If the bosons organize into a crystal, then their 
collective motions will be described by a different wave 
equation, Navier equation: dfu^ — Tl''^ djdkU^ = 0. 

Due to the particle- wave duality, our understanding of 
Euler equation and Navier equation leads to an under- 
standing of phonons in liquids and solids. We see that 
it is the organization of bosons that becomes the micro- 
scopic origin of phonons. Different organizations lead to 
different wave equations and different types of phonons. 

From this point of view, to find a microscopic origin 
of light and gravity is to find a particular organization of 
bosons (or qubits),"'^ such that the collective motions of 



such organized bosons (qubits) are described by Maxwell 
equation and Einstein equation (instead of Euler equa- 
tion or Navier equation). We like to point out that, in 
such an emergence approach, we do not look for the build- 
ing blocks of everything. The building block is known: 
it is a qubit. Even the detail form of the qubit model 
and the precise values of various coupling constants are 
not important. The important issue is how the bosons 
(qubits) are organized in the ground state. Different or- 
ganizations lead to different collective excitations, or in 
other words, different sets of "elementary" particles. 

It was shown that if bosons (qubits) organize into 
network of strings and if the string-nets form a quan- 
tum liquid (called string- net condensed state), then pho- 
tons, electrons and quarks can emerge naturally as col- 
lective motions of such organized bosons (qubits). 
In this paper, we will try to find an organization of 
bosons (qubits) such that the collective motions of bosons 
(qubits) also lead to gravitons. 

In contrast, in the traditional approach, we choose a 
particular field for each kind of elementary particle. The 
particular properties of elementary particles are encoded 
in our particular choice of the field. For example, to get a 
fermion (such as electron), we choose an anti-commuting 
Grassmann field, while to get a photon, we choose a vec- 
tor field. 

We like to point out that this paper uses a par- 
ticular emergence approach: we try to obtain every- 
thing from a local qubit model Another emergence ap- 
proach was developed in supcrstring thcory^*^ in last 10 
years, as demonstrated by the duality relations among 
various superstring models and matrix models. The 
anti-dc Sitter-space/conformal-field-theory duality even 
shows how space-time and gravity emerge from a gauge 
theory. It would be interesting find out if the grav- 
ity from the superstring theory can produce helicity ±2 
modes (gravitons) as the only gapless excitation. It is 
also interesting to find a relation between all those ideas 
of emergence. 



II. THE STRATEGY OF OUR APPROACH 
A. The rule of game 



• A practical definition of quantum gravity. 

There are many different approaches to quantum 
gravity^'^^^^ based on different principles. Some ap- 
proaches, such as loop quantum gravity,^^ stress the 
gauge structure from the diffeomorphism of the space- 
time and try to build a quantum theory where the 
whole space-time is emergent. (This is called background 
independence.^^) Other approaches, such as superstring 
theory,^"'^^ stress the renormalizability of the theory and 
try to construct a renormalizable theory that contain 
gravitons. In this paper, we will stress different things. 
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namely finiteness and locality, and follow a different rule 
of game. 

Our rule of game is encoded in the following working 
definition of quantum gravity. Quantum gravity is 

(a) A quantum theory described by a Hamiltonian in a 
Hilbcrt space. 

(b) Its Hilbert space has a finite dimension. 

(c) The Hamiltonian is a sum of local operators. 

(d) A single gapless helicity ±2 mode is the only type of 
low energy excitations. 

(e) The helicity ±2 excitations have a linear dispersion. 

(f) The gravitons^^ interact in the way consistent with 
experimental observations. 

We like to remark that different people may have different 
understandings/definitions of quantum gravity. We hope 
the above definition that will be used in this paper will 
help to avoid some possible confusions. 

Is the above definition of quantum gravity a proper def- 
inition? The first problem that we face when we consider 
quantum gravity is that "is quantum gravity a quantum 
theory?"*^ The answer to the above question appears to 
be NO, since Einstein theory tells us that time is dynam- 
ical while quantum mechanics assume an absolute time. 

The second problem is that "is quantum gravity a 
theory based real and complex numbers?" As mention 
above, the consideration of quantum theory and Einstein 
gravity suggests that the notion of manifold is not a phys- 
ical reality. Similarly, one may wonder, maybe even real 
numbers and complex numbers cannot exist as a physical 
reality. In this case, we wonder: should we avoid using 
real and complex numbers when we formulate a theory 
of quantum gravity? 

The above two considerations may serve as a long term 
goal. But they are not so useful as a guide for research. 
This is why we proposed a working definition of quantum 
gravity to guide our research at this stage. When we 
understand quantum gravity better, we expect that our 
working definition will be modified. 

Let us discuss each item in our working definition in 
more detail. Since the time is absolute in quantum me- 
chanics, the condition (a) implies that the time is not 
emergent. So the quantum gravity considered here is 
not background independent. In order for the time to 
be emergent, one has to go beyond quantum mechanics, 
that we will not do in this paper. 

The condition (b) implies that the quantum gravity 
considered here has a finite cut-off. So the renormaliz- 
ability is not an issue. A condensed matter system al- 
ways has a finite cut-off and can only simulates a system 
with a finite cut-off. So the conditions (a-b) makes it 
possible to use condensed matter system to simulate the 
quantum gravity (as defined here). Also, from a mathe- 
matical point of view, only theories with finite cut-off are 
really well defined. 

The condition (c) is a locality condition. It implies two 
additional things: (1) the total Hilbert space is a direct 
product of local Hilbert spaces Htot = ®iHi. (2) the 
local operators are defined as operators that act within 



each local Hilbert space Hi or finite products of local 
operators. The conditions (a-c) actually define a local 
bosonic model (or a local qubit model). ^ Certainly, any 
quantum spin models satisfy (a-c). It is the condition 
(d-f) that makes a theory to look like gravity. 

But why do we require gapless gravitons to be the only 
low energy excitations in the condition (d)? To answer 
this question, let us consider the theory of electromag- 
nctism where gapless photons (excitations with helicity 
±1) are the only low energy excitations. If a system also 
has helicity gapless excitations, then we are no longer 
sure if the gapless helicity ±1 excitations are photons. 
This is because phonons in solid have helicity and ±1. 
Thus a theory with gapless helicity and helicity ±1 ex- 
citations may not be a theory of photons, it may be a 
theory of phonons. But a theory with helicity ±1 excita- 
tions as the only low energy excitations must be a theory 
of electromagnetism. So here we require our theory to 
have helicity ±2 excitations as the only low energy exci- 
tations, to make sure the theory is a theory of quantum 
gravity (at least at linearized level). 

The condition (d) is a very important condition. It is 
very easy to construct a quantum model that contain he- 
licity ±2 gapless excitations, such as the theory described 
by the following Lagrangian for symmetric tensor field 

2 2 
Such a theory also contain helicity and ±1 gapless ex- 
citations and is certainly not a theory of gravity. So here 
wc impose the condition (d) rule out the above example. 
It is highly non-trivial to construct a quantum model 
that contains helicity ±2 excitations as the only gapless 
excitations. Many theories of quantum gravity fail this 
test. 

We would like to mention that, according to our defini- 
tion, only when helicity ±2 modes are the only low lying 
excitations, can the theory be a theory of quantum grav- 
ity. Such a condition may be too strict. We may want to 
relax the condition (d) to condition (d'): the helicity ±1 
and ±2 excitations are the only low energy excitations. 
In this case the helicity ±1 excitations are photons and 
the helicity ±2 excitations are gravitons, which reflects 
the situation in our universe. 

The emergent gravitons from local qubit models natu- 
rally interact with each other. However, they in general 
interact in a different way from that described by the 
higher order non-linear terms in Einstein gravity. Since 
those higher order terms in Einstein gravity are irrele- 
vant at low energies and not universal when viewed from 
the perspective of local qubit models, therefore it may 
be possible to generate those higher order terms by fine 
tuning the lattice model [such as modifying the Hamilto- 
nian {Hj and Hg), the constraints (Hu), as well as the 
Berry's phase term in eqn. (61)]. So it may be possible 
that local qubit models can generate proper non-linear 
terms to satisfy (f). 

To summarize, the goal of our approach to quantum 
gravity is to construct an Hamiltonian operator that act 
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on a Hilbcrt space. We require that the Hilbert space 
is a directly product of local Hilbert spaces which each 
has a finite dimension: H = ^iHi. Respect to such a 
locality structure of the Hilbcrt space, we require that 
the Hamiltonian operator is a sum of local operators. 
(The local operators are defined as operators that act 
within each local Hilbert space Hi or finite products of 
local operators.) In order for the Hamiltonian operator 
to describe a quantum theory of gravity, we also require 
that all the low energy excitations above the ground state 
of the Hamiltonian are described by a single mode of 
helicity ±2 excitations (gravitons). Last, those gravitons 
should interact in way that is required by the equivalence 
principle of Einstein gravity. 



B. A brief outline of our approach 



• How to design a qubit model whose only gapless ex- 
citations are described by a single helicity ±2 mode. 

We first start with a field theory of a symmetric tensor 
hij described by the following phase space Lagrangian: 

- gidkhijdkhij - g2dihijdkhkj - gzdih^jd-jhkk- 

We then put the theory on the lattice to have a finite cut- 
off. The key issue is that, as a quantum theory, does the 
lattice model has helicity ±2 modes as the only gapless 
excitation? 

It turns out that it is very hard to have helicity ±2 
modes as the only gapless excitation. In general, one 
either has all helicity ±2, ±1, and modes as gapless 
excitations or have no gapless excitation at all. In this 
paper, we show that if we compactify and descretize hij 
and its canonical conjugate tt'-' , we can actually have 
a lattice model that has helicity ±2 modes as the only 
gapless excitation. Such helicity ±2 modes correspond 
to the emergent gravitons from the lattice model. We 
also find that the low energy effective theory for such a 
lattice model has an emergent linearized diffeomorphism 
gauge symmetry /i^i. -> h^^ + + ^vf^n which is an 
implied consequence of having helicity ±2 modes as the 
only gapless excitation. 



C. Some previous approaches 

• There are many previous approaches to quantum 
gravity. But only a few of them produce a local 
quantum Hamiltonian (which satisfies the condi- 
tions (a - c)). Among those that produce a local 
quantum Hamiltonian, non of them were shown to 
contain a single helicity ±2 mode as the only gap- 
less excitation. 



In this section, we will discuss some previous ap- 
proaches in terms of our practical definition of quantum 
gravity. Superstring theory^°'^-^ satisfies the conditions 
(a), (e) and (f), but in general not (d) due to the presence 
of dilatons (massless scaler particles). The superstring 
theory (or more precisely, the superstring field theory) 
also does not satisfy the condition (b) since the cut-off is 
not explicitly implemented. 

Many approaches to quantum gravity are based on 
quantizing the classical Einstein action. The lattice grav- 
ity approach based Regge calculus on 4D space-time 
lattice^*'^^ belong to this type of approach where a fi- 
nite space-time cut-off is introduced. Such an approach 
may not satisfy the conditions (a-c) since it may not pro- 
duce a local Hamiltonian and local Hilbert space. Such 
an approach eventually failed due to the doubling phe- 
nomenon of the gravitational modes, ie it does not satis- 
fies the condition (d). 

In Ref. 33-35, quantum gravity is studied in terms of 
lattice Hamiltonians with a continuous time (just like this 
paper). However, those Hamiltonians obtained by quan- 
tizing the classical Einstein action are non-local {ie they 
do not satisfy the condition (b,c)). Also, the low energy 
excitations of the those lattice Hamiltonians were not 
discussed. We do not know if those lattice Hamiltonians 
have gapless excitation or not. As a result, we do not 
know if those lattice Hamiltonians give rise to emergent 
diffeomorphism gauge symmetry at low energies. The 
spin network'^'^ or the quantum computing^^ approach to 
quantum gravity satisfies the condition (a,b) or (a-c). 
But again, the properties (d-f) remain to be shown. 

The induced gravity from superfluid '^He discussed in 
Ref. 31 does not satisfy the condition (d) due to the pres- 
ence of gapless superfluid mode. In Ref. 32, it is proposed 
that gravitons may emerge as edge excitations of a quan- 
tum Hall state in 4 spatial dimensions. Again the condi- 
tion (d) is not satisfied due to the presence of helicity ±5* 
modes where S is unbounded. In fact, there are infinite 
many gapless modes with various helicities. 

In Ref. 36, a very interesting spin model is constructed. 
The model satisfies the condition (a - c) and 
(d" ) the gapless helicity and helicity ±2 modes are the 
only low energy excitations; 

(e" ) the helicity ±2 modes have a quadratic dispersion. 
The model is interesting since its ground state is a new 
state of matter - an algebraic spin liquid. Such a state 
is beyond the Landau's symmetry breaking description. 
Despite the emergence of gapless helicity ±2 modes, due 
to the properties (d") and (e"), the low energy proper- 
ties of the model are not very close to those of Einstein 
gravity 

In this paper, we will try to fix the two problems and 
try to construct a qubit model that satisfies the condi- 
tion (a-d), and hopeful (f) through fine tuning. We only 
partially achieve our goal. 

We studied two quantum spin models (or qubit mod- 
els). In the first quantum spin model (called the L-type 
model), the helicity ±2 gapless excitations are reliably 



7 



shown to appear as the only type of low energy exci- 
tations (ze the conditions (a-d) are satisfied). Within 
a perturbative calculation, the dispersion of the gaplcss 
helicity ±2 is found to be efc oc |fep. So the condition (e) 
is not satisfied. The appearance of the gapless helicity 
±2 modes suggests that the ground state of the quantum 
spin model is a new state of matter. 

In the second model (called the N-typc model) the 
collective modes are strongly interacting and there is no 
reliable approach to understand its low energy dynam- 
ics. Using a spin- wave/quantum- freeze approach (which 
is shown to reproduce the correct emergent C/(l) gauge 
theory in a quantum rotor model) , we argue that the sec- 
ond model may contain helicity ±2 gapless excitations as 
the only type of low energy excitations with a linear dis- 
persion u! cc k. More reliably numerical calculations are 
needed to confirm that the type-N model really has u (x k 
helicity ±2 excitations as the only low energy excitations. 

We believe that the gaplessness of the helicity ±2 exci- 
tations in both models is topologically robust: any trans- 
lation invariant perturbations cannot generate a gap for 
those helicity ±2 excitations. 



III. REVIEW OF U(l) GAUGE THEORY 

In this paper, we are going to construct a qubit model 
with emergent helicity ±2 gapless modes from a theory of 
symmetric tensor. We will show that by imposing some 
"constraints" through certain spin interaction terms, we 
can obtain a theory where the helicity ±2 excitations are 
the only low energy excitations. Since the constraints 
that we will impose are similar to the Gauss constraint 
in C/(l) gauge theory, here we will first give a brief review 
of quantum U{1) gauge theory, to introduce the physical 
ideas behind our construction in a more familiar setting. 

To obtain a.U{l) gauge theory, we may start with con- 
tinuum quantum field theory of vector fields a.i and f . 
The phase space Lagrangian has a form 

L = -S'daa^ - ^J£'£' - ^gidiajd^aj - ^g-idiUidjaj (1) 

From the equation of motion 

Sgfli = Jgid^ai + Jg2didjaj (2) 

we find the low energy excitations to be helicity ±1 and 
helicity modes with linear dispersions uj fx k. Despite 
the presence of gapless helicity ±1 modes, the above the- 
ory is not a theory of electromagnetism due to the pres- 
ence of gapless helicity mode. So the key to obtain a 
theory of electromagnetism is to gap the helicity mode. 



A. Removing helicity mode though Gauss 
constraint 



• In the standard field theory approach, we obtain 
photons by starting with a vector field theory, and 
then removing the helicity mode by imposing a 
U{1) gauge symmetry and the resulting Gauss con- 
straint. 

• The Gauss constraint makes the Hilbert space non- 
local (ie violate the condition (c)). 

The standard way to remove the helicity mode at low 
energies is to impose Gauss constraint 

d.,£' = (3) 

In the constraint system, a.i becomes a many-to-one label 
of the physical states. Different vector fields related by 
an local transformation 

ai-^a, + dif (4) 

actually label the same state. This is the well known 
gauge transformation which is generated by the Gauss 
constraint. Physical quantities should be invariant under 
such transformation so that the same state always has the 
same values of physical quantities. 

In particular, the Lagrangian for Oi field should be 
gauge invariant. An easy way to do this is to write down 
the Lagrangian in terms of gauge invariant fields. The 
gauge invariant fields are the magnetic field 

F = e'^'^dja^ (5) 

and the electric field The Lagrangian is then 

C= -£'doa.,-'H[a^,£') (6) 

n^\[j{£r + g{B'f] (7) 

The resulting equation of motion contains only two 
transverse modes corresponding to the helicity ±1 ex- 
citations. The Gauss constraint removes the helicity 
excitations from the low energy spectrum. 

However, removing helicity excitations through the 
Gauss constraint has one problem. The Hilbert space of 
the physical states is modified by the constraint, since 
only f * that satisfies the Gauss constraint are physical 
and only gauge inequivalent ai correspond to different 
physical states. The new Hilbert space can no longer 
be written as a direct product of local Hilbert spaces, 
since an arbitrary local change of in general violate 
the constraint. Thus the resulting quantum system is no 
longer a local qubit system. 

B. Try to gap helicity mode though energy 
penalty 

• Using an energy penalty to impose the Gauss con- 
straint fails, since it fails to gap the helicity mode. 
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To fix this problem, here we choose not to impose 
the constraint and not to change the Hilbert space. We 
choose instead to include a term of form U{di£^)^ in the 
Hamiltonian. The resulting theory is described by 

c = -s'doa, - H(a„r) - ic/(a,r)' (8) 

We hope that new term will suppress the fluctuations 
that violate the Gauss constraint and will gap the helic- 
ity excitations. However, from the equation of motion 
obtained from (8) 

d'^a, - g{J - Ud^){d^S,j - d,dj)a,, (9) 

we find that the extra i?7(9i£')^ term cannot gap the 
helicity mode. The helicity mode remains to have 
zero velocity. 

C. Another failed attempt 

• Imposing Gauss constraint through an energy 
penalty fails to gap helicity mode even on lattice. 

Next, we will put the continuum theory (8) on lattice 
and examine if the U{di£^)^ term can gap the helicity 
mode on lattice. To put the theory (8) on a cubic lattice, 
we introduce aij and £ij for each link (ij) of the cubic 
lattice. Here i labels the sites of the cubic lattice and Oij 
and £ij satisfy aij ~ —aji and £ij ~ The phase 

space Lagrangian for physical degrees of freedom a^j and 
£ij is given by 

{ij) (ij) (ijkl) i 

Bijki = aij + ajk + flfci + an, 

Qi= E (10) 

j next to i 

where sum over all sites, 'Yliij O'^^^ links, and 
Sijfc; sum over all square faces of the cubic lattice. This 
phase space Lagrangian tells us that aij and £ij form 
a canonical momentum-coordinate pair. The U Q? 
term is the lattice version of the U{di£'^)^ term. The 
9^ijki term and the J^i^ij^ £ij term correspond to the 
(K*)^ term and the term in the continuum theory. 

The lattice equation of motion can be easily derived from 
the Lagrangian (10). Wc find that there are still three 
gapless modes with helicity and ±1 in the long- wave- 
length limit. The U 'Y^^ term cannot gap the helicity 
mode even on lattice. 



D. Gapping helicity modes on lattice with 
discretized £ij (or campactified aij) 



• Imposing Gauss constraint through an energy 
penalty can gap helicity mode, if we (i) put the 
theory on lattice, and (ii) compactify the vector field 

All those failed attempts reveal the difficulty of gaping 
the helicity mode without gaping helicity ±1 modes 
nor changing the Hilbert space. In the following, we will 
show that in order to gap the helicity mode within the 
same local Hilbert space, we must (A) put the theory (8) 
on lattice (which is done in eqn. (10)) and (B) discretize 
£ij (or compactify Oij). 

How to discretize £ijl From (10) we see that for a fixed 
link ij, if we view aij as the coordinate of a particle on 
a line then £ij is the momentum of the particle. To dis- 
cretize (or quantize) the momentum £ij, we simply put 
the particle on a circle instead of a line. This is achieved 
by letting aij and Oij + 27r to describe the same point 
(thus to compactify Oij). After the compactification, £ij 
is quantized as integer and the phase space Lagrangian 
(10) need to be modified to 

i = E ^ijdoaij - ^-^E ^ij + -9 E cos(Bijfei) 

(ij) (ij) (ijkl) 

-^uY^Ql (11) 

i 

in order to be consistent with the periodic condition 
Oij ^ aij + 2tt. (11) describes a rotor model with ro- 
tors on the links of a cubic lattice. It was shown that in 
the U ^ g ^ J limit, the ground state of the rotor model 
is a string-net condensed state. The low energy ex- 
citations above the string-net condensed state are shown 
to be gapless helicity ±1 modes and the helicity mode 
is gapped! ^^'^^"'^^ Basically, after quantization, the oper- 
ator has a discrete spectrum. So in the large U 
limit, the term U J^i Qi remove some fluctuations from 
the low energy spectrum. The removed fluctuations turn 
out to be the helicity mode. 

E. Low energy collective modes of the rotor model 
through spin- wave/quantum- freeze approach 

• When we treat our lattice model as a classical the- 
ory, we find some classical gapless modes have very 
weak quantum fluctuations, and they remain gap- 
less in quantum theory. While other classical gap- 
less modes have very strong quantum fluctuations, 
and they acquire a gap in quantum theory. 

The gapping of the helicity mode can also be un- 
derstood from a spin- wave approach if we incorporated a 
mechanism called quantum freeze. We will use a similar 
approach to study the emergence of helicity ±2 gapless 
mode. 

To understand the low energy dynamics of the com- 
plicated and strongly interacting rotor model (11), let 
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us treat the model as a classical model and {aij,£ij) as 
classical fields. The classical ground state is obtained by 
minimizing the Hamiltonian 

(ij) (ijkl) i 

We find that the classical ground state is given by 
{aij,£ij) — (0,0). The classical low energy collective 
modes is given by the fluctuations aij and £ij . 

One way to obtain the dynamics of the classical collec- 
tive modes is to obtain the continuum effective theory of 
lattice model (11). To obtain the continuum theory, let 
us assume the fluctuations of aij are small and expand 
(11) to the quadratic order of aij. Then we take the con- 
tinuum limit by introducing two vector fields (a, £) and 
identifying 

aij = f da;*a*, £ij = f dx"£\ 

J i J i 

Bijki = [ dS-idxa) 

J ijkl 

where /^^-^j dS is the surface integration on the square 
ijkl, and we have assumed that the lattice constant a ~ 
1. The resulting continuum effective theory is given by 
(6). The classical collective modes described by a and 
£ are the spin-wave-like fluctuations in the rotor model. 
We find that there are three gapless modes with helicity 
and ±1 at classical level. 

Now the question is that do we trust the above classical 
spin wave result? So in the following, we will study the 
quantum fiuctuations of those classical modes to check 
the self consistency of the classical spin- wave approach. 

To study the quantum fluctuations of £'' and a^, we 
note that the longitudinal mode and the transverse 
modes separate. Introduce £ = £\\ +£± and a = a\\ +a±, 
we find that the dynamics of the transverse mode is de- 
scribed by 

C± = a±do£± - ^£1 - ^dia±dia± 

At the lattice scale 6x^1, the quantum fluctuations of 
£± and a± are given by 

We see that when J <^ g the fluctuations of a is much 
less than 1. So expanding cos{Bijki) to quadratic order 
is a good approximation. But this alone does not grantee 
the validity of the spin-wave approach. In quantum the- 
ory, the compactness of aij imply that £ij is discrete. 
So we cannot treat £ij as a continuous variable, as we 
did in the classical spin-wave approach. However, in the 
J <^ g limit, we see that the quantum fluctuations of the 
£ (or £ij) is much larger than 1 which is the discreteness 
of £ij . In this case, we can indeed treat £ij as continu- 
ous variables. So the spin- wave approach is valid in the 



J ^ g limit for the transverse mode. The classical spin 
wave result can be trusted even in quantum theory. We 
conclude that the transverse mode (or the helicity ±1) 
has a linear gapless dispersion. 

The longitudinal mode is described by two scalar flelds 
(/(a;), 7r(a;)) with = dif and tt = di£^ . Its dynamics 
is determined by 

^11 ='rao/-^^(-9-2)7r-i[/^^ 

At the lattice scale, the quantum fluctuations of tt and 
/ are given by Jtt = and Sf = 00. We see that the 
fiuctuations of / are much bigger than the compactifica- 
tion size 27r and the fiuctuations of tt are much less then 
the discreteness of f ' which is 1. In this limit, the result 
from the classical spin-wave approach cannot be trusted 
once we include quantum effect. 

In fact the weak quantum fluctuations in the discrete 
variable tt and the strong quantum fiuctuations in the 
compact variable / suggest that the corresponding mode 
is gapped after the quantization. Such a phenomenon is 
called quantum freeze. 

Since tt has weak fiuctuations which is less than the 
discreteness of tt, the ground state is basically given by 
TT = 0. A low lying excitation is then given by tt = ev- 
erywhere except in a unit cell where tt = 1. Such an ex- 
citation have an energy of order U. The gapping of helic- 
ity mode is confirmed by more careful calculations. ^^'^^ 
From those calculations, we flnd that the weak fluctua- 
tions of TT lead to a constraint tt — di£'^ — and the 
strong fluctuations of / lead to a gauge transformation 
ai — > ai + dif for the low energy states. The Lagrangian 
(8) equipped with the above constraint and the gauge 
transformation becomes the Lagrangian of a U{1) gauge 
theory. 

To summarize, the emergence of photons (the helicity 
±1 excitations) from a local qubit model is purely a quan- 
tum effect which requires a discretization of the electric 
field on lattice. It also requires the ground state of the 
qubit model to have a new kind of order - string-net con- 
densation (which is ensured by the condition J <C p). In 
this case, the helicity ±1 modes have small fluctuations 
and are gapless, while the helicity mode has extremely 
strong quantum fluctuations and are gapped. This ex- 
ample gives us important hints on what are required in 
order to have emergent helicity ±2 excitations from a 
qubit model. 



IV. GENERALIZATION TO SYMMETRIC 
TENSOR FIELD 

A. A field theory of symmetric tensor 



• The phase-space Lagrangian of a symmetric-tensor 
field theory. 
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To construct a qubit model with emergent helicity ±2 
gapless modes, let us start with a continuum field theory 
with symmetric tensor fields aij and f 

where i^j = 1,2,3 are the spatial indices. The phase 
space Lagrangian of the field theory is given by 

- gidkaijdkaij - g-idiaijOkakj ~ gsdiaijdjakk- 

We see that is the canonical momentums of aij and 
satisfy 

[a,j(y),r^^{x)] = ^(<5,™<5,„ + - y) (12) 

as operators. 

The model can be put on a lattice and becomes a qubit 
model. From the resulting equation of motion, we find 
that there are six gapless linear modes. Two modes with 
helicity and four modes with helicities ±1 and ±2. So 
to obtain emergent graviton, we simply need to gap the 
two helicity modes and the two helicity ±1 modes while 
keep the helicity ±2 modes gapless. 

B. di£^-' — constraint and helicity and ±2 
quadratic modes 

• The vector constraint and the associated gauge sym- 
metry remove the helicity and ±1 modes in sym- 
metric tensor field theory. 

First, let us try to remove those modes through con- 
straints in the continuum theory. We first impose three 
constraints 

d.,£'^ = (13) 

which are called the vector constraints. Obviously these 
three constraints commute with each other. From the 
constraint, we can construct a unitary operator 

that generates gauge transformations on aij : 

a,j WoijW'^ = a,j + d.Jj +djft (15) 

The gauge invariant field is a symmetric tensor field 

= Ri^ = e™'^-e^'"a™aza„fc (16) 

The dynamics of the constraint system is describe by 
a gauge invariant Lagrangian. To the lowest order, such 
Lagrangian has a form 

C = £'^doa,^ - a{£'^f - P{£''f - 7{R'^f - A(i?")^17) 



with q;, /3,7, A > 0. After solving the equations motion, 
one finds that the constraint system contains three gap- 
less modes with quadratic dispersion w oc fc^. Two modes 
have helicities ±2 and one has helicity 0. Thus the vector 
constraints (13) remove the two helicity ±1 modes and 
one helicity mode. 

In fact, this constraint model defined through (13) and 
(17) is the continuum limit of the lattice model studied in 
Ref. 36. Due to the quadratic dispersion and the helicity 
gapless mode, the lattice model studied in Ref. 36 does 
not reproduce Einstein gravity at low energies. However, 
the lattice model does give rise to a new quantum liquid 
of qubits (or quantum spins). Just like any other ordered 
state, the new quantum liquid also contain topological 
defects. As in U(l) gauge theory, we can allow the three 
constraints (13) to be violated at one point 

d^£'^x) = q^6{x-xo) (18) 

and create a defect. Such a defect is called "electric" 
charge. It is clear that the "electric" charge g' is a vector. 
We also note that there are three identities for i?*^ 

^^R'^ = (19) 

in the continuum. Those identities can be violated on 
lattice. The violation of these identities will also create a 
local defect called the monopole. Again the correspond- 
ing "magnetic" charge is a vector. 

In Einstein gravity, the mass of matter is coupled to 
gravity via a modification of constraint. In other words, 
the way that mass generates gravity is similar to the way 
that electric or magnetic charges generate electromag- 
netism. From the above discussion we see that the model 
defined through (13) and (17) has no scaler charges (elec- 
tric or magnetic) that may correspond to mass. This is 
a more fundamental reason why the model can not be a 
theory of quantum gravity. 

C. Additional R" — constraint to remove the 
second helicity mode 

• The additional scaler constraint and the associated 
gauge symmetry remove the last helicity mode in 
symmetric tensor field theory. 

• In the constrained model, the gapless helicity ±2 
modes have a ui ^ dispersion. 

• The scaler constraint can be violated at an isolated 
point, which corresponds to a mass. 

To obtain a theory with emergent gravitons and point 
scalar mass, we would like to remove the second helicity 
mode by introducing another constraint. This constraint 
must be a scalar constraint. There are only two scalar 
physical constraints = and i?** = that can be 
implemented in the theory. 
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If wc use the constraint 

= 



(20) 



then the i?*^ is no longer gauge invariant under the gauge 
transformation generated by the new constraint. The 
lowest order gauge invariant tensor field are third order 
derivatives of aij . The constraint does remove the helicity 
modes. The helicity ±2 modes have a cubic dispersion 
relation uj (x . Furthermore, we find that the "elec- 
tric charge" cannot even be point-like. Only electric line 
charge is allowed in this case. We see that the resulting 
theory is very different from Einstein gravity. So we will 
not go further in this direction, although the resulting 
qubit system is a very interesting condensed matter sys- 
tem which gives rise to a new quantum state of matter. 
The right constraint that we should use is 



W = 



In terms of a^ , this constraint is 



[Si^d"^ - didj)aij = 



(21) 



(22) 



Violation of this constraint at one point corresponds to a 
scalar charge which can be interpreted as a point mass. 

But this constraint is very unusual, it is a local con- 
straint on aij itself, so it generates a gauge transforma- 
tion on 



(23) 



Under such a gauge transformation, a linear gauge in- 
variant field is given by 



(24) 



Because aij is the conjugate variable of , wc should 
carefully check the commutating relation of these four 
constraints. Wc first transform aij, into Fourier space 
and it is easy to check 

K„(fc),r^(fe')] = '-[s^mSjn + s,jj„^]i27Tfsik - k') 

(25) 

then we only need to check 

[fcj£ -'(fe), {dnmk fc„fc„j)a,j,n(fc )] 

= '-i2n)%i5„^k" KXrMnAn + 5.„5,,„]<5(fc - k') 



i{2TTykj{S,jk" - klk'j)S{k - k') 




(26) 



It turns out that these four constraints commute with 
each other. 

To obtain the dynamical properties of the symmetric 
tensor field theory with the constraint (13) and (21), we 
need to find a gauge invariant Hamiltonian. One way 



to construct a gauge invariant Hamiltonian is to express 
the Hamiltonian in terms of gauge invariant fields Cj in 
eqn. (24) and i?'-'. This way we obtain the following 
Hamiltonian density 



2 ■' •' 2 
or the Lagrangian density 

£ ^ £'Woa^J - ^qC} - R^' . 

We will call such a system L-type model. 

The equation of motion for this Hamiltonian is 



(27) 



(28) 



£^^ = -.ge"™^e'='^e""''='e^''"'c>„,9zc>„,.9pa„a-' (29) 

If we use a plane wave solution along the third direction, 
with fci = 0, fc2 = 0, fcs = fc, it is easy to find 





= gJk' - 










£'' 








(30) 

We see that there are only two modes of gapless excita- 
tions, — g^Yid £^^, with Ld (X k^ dispersions. The 
fluctuations of form £^^ + £22 j^^g 2cro frequency and 
represent the pure gauge fluctuations of (23). So only 
the two helicity ±2 modes are physical and the second 
helicity mode is removed. 

We can further calculate the force between two scalar 
defects (the masses) obtained by violating the constraint 
(22): 



R"{x) = mid{x — Xi) + m2S{x — x-z) 



(31) 



We find the force is repulsive between two masses with 
the same sign and decay as \xi — X2\^'^- 

To summarize, we constructed a L-type model using 
symmetric fields that satisfy the commutation relation 
(12) and the constraint (13) and (21). The gauge invari- 
ant Hamiltonian is given by (32). The only low energy 
excitations of the L-type model are gapless helicity ±2 
modes with a cubic dispersion. Although the helicity ±2 
excitations are the only low energy excitations, the low 
energy properties of the L-type model are very different 
from those of Einstein gravity where the helicity ±2 ex- 
citations have a linear dispersion. 



D. Nonlocal gauge invariant Hamiltonian 
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• In the constrained model, the gapless helicity ±2 
modes can have a linear uj ^ k dispersion, if the 
field theory Lagrangian is gauge invariant only up 
to a surface term. 

How to construct a model that has a linearly dispersing 
graviton and 1 /r^ attractive force between point masses 
with the same sign? One way to achieve this is to choose 
a nonlocal gauge invariant Hamiltonian (i.e. the Hamil- 
tonian density is gauge invariant up to a total derivative). 
Let us consider the following Hamiltonian density 



or the Lagrangian density 



The quadratic term 



(32) 



(33) 



(34) 



is gauge invariant under (23), provided that £^^ satisfy 
the constraint (13) di£'^^ — 0. From the identity (19), we 
also find that OijK'^ is invariant under the gauge trans- 
formation (15) up to a total derivative. Thus the total 
Hamiltonian is invariant under the gauge transformation 
(23) and (15). 

The corresponding equations of motion have the form 





= -gJk'(^£''- 






= -gJk' (^£" - 


¥ 




= .gJfc2£i2 






= £23 ^ £33 ^ Q 





(35) 



From the solution of the above equations of motion, we 
find that the only gapless modes are the two helicity ±2 
modes with a linear dispersion relation. Those modes 
can be identified as gravitons after quantization. 

We can create two point defects through (31). We find 
that the force between the two masses is proportional 
to TOi77i2/r2. Two masses with the same sign attract 
and two masses with opposite signs repel. The quantum 
gravity theory defined through the commutation relation 
(12), the constraint (13) and (21), and non-local gauge 
invariant Hamiltonian (32) is called N-type model. It is 
the N-type model that represent a quantum theory of 
Einstein gravity (at the linear level). 

We would like to point out that when compared to 
Einstein theory of gravity, Oy- can be interpreted as the 
fiuctuations of the spatial part of a metric tensor g^^, 
around fiat space: ~ gij — Sij. is related to the 
three dimensional Ricci tensor 7^*^ : 

1 



If we introduce ooo and a^t as Lagrangian multipliers 
to impose the vector and the scaler constraints, wc can 
rewrite our model (27) as the following Lagrangian 



C = £^^doaij — ^ 



1 2 



+2aQidj£^^ + 000(92 - didjOij) 



(37) 



(36) 



After integrating out £'■', we find that this action is ex- 
actly the linearized Einstein action around a flat space- 
time. The associate gauge transformations in space are 
also enlarged to gauge transformations in space-time. 



THE LATTICE MODELS AND IMPOSING 
CONSTRAINTS THROUGH ENERGY 
PENALTIES 



• Put the symmetric-tensor field theory on lattice and 
construct qubit models. 

So far, we have discussed continuum models that has 
gravitons as the only gapless excitations. The gapless 
helicity and ±1 modes are removed by imposing the 
constraints (13) and (22). However, the two constraints 
changes the Hilbert space of the model. The new Hilbert 
space docs not have a local form Ti = ®„H„ where UnS 
are the local Hilbert spaces. Thus the constrained models 
are not local qubit models. In this section, we are going 
to fix this problem. We will show that the vector and 
the scalar constraints can be realized in lattice models 
through certain energy penalty terms without changing 
the Hilbert space. 

First let us put the symmetric tensor field theory on a 
cubic lattice. Ref. 36 propose a nice way to do so. Here 
we will follow that convention to put Oxx, Oyy, Ozz, £^^, 
£yy, and £^^ on the vertices and put Qxy, ctyz, azx, £^^, 
£y^, and £^^ on the square faces of the cubic lattice. For 
example 

axx{x) Oxxii), 

azz{x) -> azz[i), 

axy{x) axy{i+- + -), 

, , , . y z 

ayz{x) ayz[t + - + -), 

Ozxix) azx{i +2^2'-'' 

where i is the integral vector that represents the position 
of a vertex of the cubic lattice. 



A. Putting the constraints on lattice 

• The vector constraint and the scaler constraint are 
put on lattice via an energy penalty term. 
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On the cubic lattice, one of vector constraints in 
eqn. (13), diE"^^ , can be written in the foUowing form 

Q{i, i + x) = + x)- S'-'^i) 

+£:-(i + | + f)_f-(i + |_£) 

The other two constraints become Q{i, i+y) and Q{i, i + 
z) which are obtained from the above expression by cy- 
chng xyz to yzx and zxy. We note that the constraint 
field diS''^ becomes a quantities on the links in the j- 
direction. The scalar constraint (22) becomes 

Vii) = X! [abb{i + a) + abb{i - a) - 2abb{i)] 

a,b—x,y,z 

- ^ [Uaaii + a) + aaaii - a) - 2aaa{i)] 

r / . o, b . . , a b . 

-2 X 2^ [a„b(^ + - + -)- aa^,{^ - 2 + 2^ 

ab—xy.yz^zx 



,. a b ,. a b 

Let us introduce, L*-' and 9ij. For example 

i^^(i) = (-)^£:^^(i) 

i^^i + f + f) = -(-r^^ni + f + f) 



and 



Oxx{i) = {-Yaxxii) 



U{i + f + f ) = -2(-)'a.-;(i + ? + ?) 



2 2- 



L"'' and 6ab satisfy the following commutation relation 
between angle and angular momentum: 

[L^^ii), Oabii)] = -i, others = 0. 

Using and dab, we can rewrite the above discretized 
constraints as 

Qii,i + x) =i-""(i + a;) + L""(i) 

^ 2 2' ^ 2 2' 

^ 2 2' ^ 2 2' 



and 



77(i) 



^ [6lt6(i + a) + 6ibb(i-a) + 2e'fcf,(i)] 

^ [6laa(* + a) + 6'aa(i-a)+26'aa(i)] 



a— a:,y ,. 



Ea b ^ , a b 

[^ab(* + 2 + 2) + ^"''(* - 2 + 2^ 



-^Qab{t+^-^)+eab{t-^-^ 



One can check that (5(i, i + a) and 7](j) all commute with 
each others. 

Just as discussed in section III, we would like to impose 
the r](%) = and Q(i, i + a) =0 constraints by including 
the term 



in the Hamiltonian without changing the Hilbert space. 

However, using H'jj to impose the r](%) = constraint 
docs not work. This is because 77 has a continuous spec- 
trum. No matter how large is [/21 the excitations that 
violate the r\(%) = constraint do not have an energy 
gap. So the low energy excitations do not have to satisfy 
the r]{i) = constraint. Similarly, if Qab are not compact 
and L°-^ are not discrete, the Q{i,i + a) operators will 
also have continuous spectra. The low energy excitations 
do not have to satisfy the Q{i, i + a) = constraint. 



B. Compactification 



• To use the energy penalty to impose the vector and 
scaler constraints, we need to compactify and de- 



scretize the lattice fields a^x, o.yy, a^z, axy, Oyz, -^zx- 



a. 



In order to use the energy penalty H'^ to impose the 
r{(%) = and Q{i,i + a) = constraint, we need to 
make the spectra of ri{i) and Q{i, i + a) discrete. In this 
case, for large ?7i,2, all excitations below an energy gap 
of order Ui^2 will satisfy the ri{i) = and Q{i, i + a) = 
constraint. To discretize the spectrum of the constraint 
operators ry(i) and Q{i, i + a), wc need to discretize both 
L'^^ and 6ab or compactify both L"'' and Oab- 

We compactify 9ab by imposing a periodic condition 
^ Oab + 27r. Similarly, we compactify L"'^ by impos- 
ing a periodic condition L°-^ ^ L"-^ -f no where uq is 
an integer. 6ab and are no longer physical operators 
after the compactification. Only Wl^ = (P'^^L'^* /'"■g ^nd 
Wg^ = e*^"'' and their products are physical operators. 
For a fixed ab and i, W^^{t) and VF^''(i) satisfy the al- 
gebra 



Wf{i)Wt(i) 



''''''Wf{i)Wl\i) (38) 



Such an algebra has only one uq dimensional represen- 
tation. This no dimensional representation becomes our 
local Hilbert space Hi^ab- The total Hilbert space is given 
by Ti = ®i.abT~Li.ab after the compactification. In other 
words, there arc states on each vertex and uq states 
on each square face of the cubic lattice. The total di- 
mension of the Hilbert space is finite for a finite cubic 
lattice. Wc note that in the uq 00 limit, we recover 
the uncompactified case. 

The constraint operators r](i) and Q(i,i + a) are not 
allowed for the campactified model since they are not 
products of Wl'' and Wg"''. This can be fixed easily. In 
the campactified model wc replace ?7^(i) and Q^{i,+a) 
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hy l — cos[ri{i)] and nQ{l — cos[27TQ{i,i + a)/nG]) respec- 
tively. So, in the compactified model, the terms in the 
Hamiltonian that impose the constraints have the follow- 
ing form 



(39) 



+ nGUiY^ {l-cos[2TTQ{i,i + a)/nG]} 

i a—x,y,z 

The operators cos[77(i)] and cos[2TTQ{i,j)/nG] have dis- 
crete eigenvalues. We will also show below that those 
operators all commute with each other. So the Hamilto- 
nian Hu has a finite energy gap. The low energy states 
below the gap (i.e. zero energy states) all satisfy 



_ g2i7r(5(i, i+o)/riG 



1. 



(40) 




We hope that Hu in the compactified model can gap all 
the helicity and ±1 modes. 



C. Low energy Hilbert space 



• The structure of the low energy Hilbert space of the 
lattice model that satisfies the vector and scaler con- 
straints Hjj = 0. 

In large Ui and IJ2 limit, the low energy Hilbert space 
satisfy the scalar and vector constraints 



S{i) 



V{i,i- 



cx:p[iri{i)] = 1 
cxp[27riQ(i, i - 



a)/nG] 



1 



(41) 
(42) 



where a = x, y, z. The scalar constraint has the following 
form 



(43) 



The structure of the integer coefficients c?^^ can be seen 
more clearly in Fig. 1. Fig. 1 also shows the action of 
S{i). The vector constraint has the following form 



V{i, i + x) = exp(i 



(44) 



The structure of the integer coefficients d'}\^^ ^ and the 
action of V{i, i + x) are plotted in Fig. 2. 

From the Figs. 1 and 2, one can see that 
Ci,r'^j'j+x,r = 0' where sum over ah the ver- 
tices and the square faces of the cubic lattice. So the 
scalar and vector constraints commute with each other, 
[S{i),V{j,j + a)] = a = x,y,z. We can construct 
states that satisfy all the constraints one by one. 

The states of the spins on vertices and squares can be 
labeled as 



FIG. 1: (Color online) The integer coefficients c°'''{i, r) in (43). 
A short solid line in the a direction marks a non-zero c°'°'{i,j). 
The thicker lines correspond to c°''^{i,j) — 4 and thinner lines 
c'^°'{i,j) = 1. The filled circles mark the non-zero c^''{i,i + 
^ + |) etc whose values are equal to —1. Red represents 
positive integers and blue negative integers. The action of 
the scaler constraint operator S{i) in eqn. (41) changes L"'''{r) 
mod no- The figure also shows those changes. The short solid 
lines mark the non-zero changes of L"" (i) where thicker lines 
correspond to a change of 4 and thinner lines a change of 1. 
The filled circles mark the non-zero changes of L^''(i + j + ^) 
etc . The changes are equal to —1. 



l+Z 







• 






• 

1 



l+X 



FIG. 2: (Color online) The integer coefficients d'^''{i, i + x,r) 
in eqn. (44). The short solid lines mark the non-zero d°'°'{i, 
x,j) which are equal to 1. The filled circles mark the non- 
zero d^^{i, i + x,i + ^ + ^) etc , which are also equal to 1. 
The action of the vector constraint operator V{i, i + x) in 
(42) changes nG0"'^(r)/2Tv mod no- The figure shows those 
changes. The short solid lines mark the non-zero changes of 
nG0°'°'{i)/2-K. The filled circles mark the non-zero changes of 
nod'^^ii + J + f )/27r etc . All the changes are equal to 1. 



0,l,2---nG-l 



(45) 
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which satisfy |n) = |n + uq) and 

Wl\i)\n)i^ab = e'^^^'''^\n)i,ab 
Wt{.i)\n)i^ab = \n+l)i,ab (46) 

it is easy to show {Wg''{i)Y° \n)i^ab = \n)i^ab- 

The simplest state that satisfy the vector constraints 
(42) is the state |0) = <^i ab |0)i,ab- However, this 
state does not satisfy scalar constraint (41). Note that 
S{i)\Q) is a new state still satisfies the vector constraints 
(42). Since {S{i)f |0) satisfies (42) and 5'(i)"'^|0) = |0), 
therefore the state |6) = (g)- {tZ%S{iY) |0). satisfies 
both the vector constraints (41) and scalar constraint 
(42). Similarly, starting from any state that satisfies 
the vector constraints (42), we can apply the operator 
®i (X]n=o '5'(*)") to obtain a state that satisfies all the 
constraints. 

Now we only need to construct the states that 
satisfy the vector constraints (42). All these 
states can be generated by acting the operators 
T^x. (*) , Te^y (i) , 7^,, (i) , 7^,2, (i + a;/2 + y /2) , 7^.y, (i + y /2 + 
2;/2), 7?.22;(i + 2;/2 + a;/2) on vacuum |0). Those operators 
are actually the discretized version of i?'-* , for example: 

T^yzii + f + f ) = cxp{^[20,,(^) + 2e,,{i + y) 

+2e^^{i + z) + 2e.^.,.{i + y + z)- e^y{i + 1 + 1) 

-^xy(» - 2 + f ) " + 2 + 2 ^ " " 2 + 2 ^ 

-e,.y{i + ^ + 1 + 1) - e,y{i + 2 - 1 + 1) 
-9,^{i + y + I + |) - 0,.,{i + y - I + |) 

+0y,{i + a; + I + |) + 9y,ii - a: + I + |) 

+20,.(i + | + |)]} (47) 

T^zz{i) = cxp{i[0xx{i + y) + 6'^x(i - y) + 26*2;.^ (i) 
+0yj^(i + x)+ 9yy(i -x) + 26'j^j^(i) - 9^y{i + ^ + 2 ) 

-^..(^ - f + f ) - ^^^.(^ + f - f ) - ^-^(^ - f - f )]} 

(48) 

Other components can be easily obtained by cycling 
xyz to yzx and zxy. 

Figures 3 and 4 illustrate the action of TZab- Compar- 
ing the patterns in Figs. 3 and 4 with the pattern in 
Fig. 2, we can show that the operators TZab commute 
with the vector constraint V{i,i + a). As functions of 
9ab, the operators TZab also commute with the scalar con- 
straint S{i). The operators that commute with the scalar 
and the vector constraints are called gauge invariant op- 
erators. It is easy to check TZab are the complete gauge 
invariant operators and the low energy Hilbcrt space can 
be constructed by acting these operators on |0) one or 
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FIG. 3: (Color online) The action of the operator TZzzii) in 
eqn. (48) changes L°'''{r) mod na- The short solid lines mark 
the non-zero changes of L°'°'{i). The thicker lines correspond 
to an change of 2 and thinner hnes a change of 1. The filled 
circles mark the non-zero changes of L^^{i + § + f ) etc . The 
changes are equal to —1. 
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FIG. 4: (Color online) The action of the operator TZyz{i + 
^ + §) in eqn. (47) changes L''''{r) mod na- The short solid 
lines mark the non-zero changes of L"''^{i) where the changes 
are equal to 2. The filled circles mark the non-zero changes of 
L^^{i -I- § -I- f ) etc , where the bigger circle corresponds to a 
change of 2 and smaller circles a change of ±1. Red represents 
positive integers and blue negative integers. 

more times. But one should notice not all such oper- 
ations can successfully create new states, for example, 
TZ,^{i)TZyy{i)TZzz{i)\0) = S{i)\b) = |0). 

VI. THE L-TYPE LATTICE MODEL 



• Giving non-trivial dynamics to states in the Hjj = 
subspace through Hj and Hg . 

In the last section, we have considered spin model 
whose Hilbert space is defined through the algebra (38) 
and whose Hamiltonian is given by Hij. We find that 
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such a system has infinite many zero-energy states and all 
other states have an energy at least of order C/i,2/"g- In 
this section, we are going to add additional terms to the 
Hamiltonian that act within the zero-energy subspace. 
The zero-energy subspace is called the constrained sub- 
space and the operators that act within the constrained 
subspace are called gauge invariant operators. The new 
term will lift the degeneracy of the zero-energy states 
and give those states a non-trivial dynamics. We find 
that the resulting low energy collective modes contain 
only helicity ±2 modes with a; ~ fc"^ dispersion. 

A. Putting the term on lattice 







i+z 


• 

-1 
• 




+y ^1 





i i+x 



FIG. 5: (Color online) The integer coefficients f^^cd,r in 
eqn. (52) (see eqn. (49)). The filled circles mark the non-zero 
fx,ad,r , which have values ±1. Red represents positive 
integers and blue negative integers. 



We have seen that TZ""^ is a gauge invariant operator 



which acts within the constrained subspace. 
add the following term 



So we can 



nog 



J2 '^^^ 



h.c. 



i.a—x.y,z 



i ,ab—xy ,y z .zx 



to the Hamiltonian. It turns out that if we expand gauge 
invariant operator 7^"'' to quadratic order in 6"^^ , we will 
obtain the gauge invariant term R^^ . Thus the above 
lattice Hamiltonian term corresponds the R^^ R^^ term in 
the continuum model. 



B. Putting the C!C! term on lattice 



-1 





i+z 


2* 


• 



i+y 



FIG. 6: (Color online) The integer coefficients ^ in 

eqn. (52) (see eqn. (50)). The short solid fines mark the non- 
zero f!^^^j, which have values ±1. The filled circles mark 
the non-zero ^ etc , which have values 2. Red represents 
positive integers and blue negative integers. 



To put the C'jCj term in the continuum model on lat- 
tice, we first put the field Cj on lattice. From the relation 
between continuum field and the lattice operator L'-* , 
we find a set of lattice operators crj that correspond to 
Cj. <7j are given by, for example 

= _L-(i + | + f)-L-(i+| + |+y) 



L^^i + ^ + l + z) (49) 



where /^^^ ^ are integer coefficients. Figs. 5, 6, and 
7 show that pattern of those coefficients. From those 
patterns, one can check that the operators commute 
with the constraint operators S{i) and V{i, i + a). 

After compactification, the corresponding physical op- 
erators are defined as 



T'a = exp 



2111 



no 



(53) 



Thus Hamiltonian term CjC* 



in the continuum model 



al ^ 2Ly^i + f + f ) + ^L'^i - I + f ) 

- [Lyy{i + z)- L'^'={i + z)- L'-'{i + z)] 

- [Lyy{i)-L^''{i)-L'-\i)] (50) 

al = [L-'{i + y)-L--{i + v)-Lyy{i + y)] 
+ [L""(2) -L""(f) -Ly^'(i)] 

- 2Ly'(i + ^ + -) - 2Ly'(i + ---) (51) 

^2 2' ^2 2' ' 

We note that crj^ have a form 

< = E rici.r^'\r) (52) 
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FIG. 7: (Color online) The integer coefficients /fcd.r in 
eqn. (52) (see eqn. (51)). The short solid lines mark the non- 
zero fz,cc,jj which have values ±1. The filled circles mark the 
non-zero /fed r , which have values —2. Red represents 
positive integers and blue negative integers. 
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become the follow lattice Hamiltonian term 



i.a—x^y^z i^ab—xy^yz^zx 



Putting everything together, we obtain the following 
lattice Hamiltonian 



a 



(54) 



Such a lattice Hamiltonian defines a L-typc lattice model. 



C. Low energy collective modes — 
spin- wave / quantum-freeze approach 

• The only gapless exactions in the L-type model are 
helicity ±2 modes with uj dispersion. 

• The low energy effective theory for the L-type model 
has an emergent diffeomorphism gauge symmetry. 

In this section, we will use spin- wave/quantum- freeze 
approach to study the low energy dynamics of the L-type 
model (54). Another approach, the equation-of- motion 
approach, will be given in appendix A. 

In the spin- wave/quantum- freeze approach, we first 
assume the quantum fluctuations of 9ab and 0°'' = 
/no to be much larger than their discreteness 
Mab ^ ^(p"'' ^ I /no in the large uq limit. (We will 
check this assumption later for self consistency.) In this 
case, we can treat 9ab and t/)"^ as classical fields and use 
the classical phase-space Lagrangian (61) to study the 
low energy collective modes. 



Classical spin waves 



Let us also assume that the quantum fluctuations of 
9ab, 4'°'^ are much smaller then 1 so that we can expand 
the phase-space Lagrangian to quadratic order in 9ab, 
(f)"-^ . This allows us to find the dispersions of collective 
modes of our bosonic model (54). There are total of six 
collective modes. We find four of them have zero fre- 
quency for all fc, and two modes have a cubic dispersion 
relation near k = (tt, 7r,7r). Near k = (tt, 7r,7r) , the dy- 
namics of the six modes are described by the following 
continuum field theory: 



£ = no 



2n 



9_ 
2 

nGU2 



(i?")2 



and C] 



(55) 



where 

e™"a„ - i<5„j>"). The continuum fields 
are given by —(—1) k" ' 2 
{-lye^^ii) for a = b, where a,b 
uum fields (jf''^{x) arc given by — (— 1)''(/)' 
a 7^ 6 and by (-1)>°''(2) for a = h. 



dab{x) 

+ % + I) for a 7^ 6 and by 



x,y,z. 

kabf^ 



The contin- 
f f + |)for 



2. The classical helicity ±2 modes 

The two modes that have a cubic dispersion relation 
are the two helicity ±2 modes. The dispersion relation 
is given by Wfc ^ \J g^J^\k'^. We note that the Ui and 
U\ terms decouple from the helicity ±2 modes, since the 
helicity ±2 modes satisfy the constraints c?i</)*^ = i?*' = 0. 
Thus the dynamics of the helicity ±2 modes does not 
depend on \Ji and IJ\. 

To see if the result Uk ^ \J g^ J^\kf can be trusted, 
we need to analyze the quantum fluctuations of Qij and 
. We find that for large nc, the quantum fluctu ations 
of the hehcity ±2 modes is of order 5(j)^^ ,59ij ^ \/l/nG 
(note that and g^ are of the same order). So the 
fluctuations of (j)^^ and 9ij satisfy I /no 6(1)'^^ <C 1, 
1/riG ^ S9ij <C 1 and the semiclassical approximation is 
valid for the helicity ±2 modes. In this case, the result 



gJ\k\^ can be trusted. 



3. The quantum freeze and the gapping of helicity and ±1 

modes 

The helicity ±1 modes and one of the helicity mode 
are described by {9i, (pi) which correspond to fluctuations 
of the following form: 9ij = di9j -\- dj9i and 4>^ = dj(f>^^. 
Their frequency = 0. For such modes, the Hamilto- 
nian only contains (p' . Thus the quantum fluctuations 
satisfles Scj)'' ^ l/?^G and S9i ^1. So the semiclassical 
approximation is not valid and the result = cannot 
be trusted. Using the similar argument used in the emer- 
gence of U{1) gauge bosons in section HID, we conclude 
that those modes are quantum frozen and are gapped. 
Since discreteness of 4>^^ is of order l/nc, the gap of the 
helicity and ±1 modes is of order Ui/na. The strong 
fluctuations S9ij = di9j + dj9i 1 lead to gauge trans- 
formations 

9,j ^ 9,, + d,9, + d,9, (56) 
and the weak fluctuations 0* <C 1 /no lead to constraints 







(57) 



The second helicity mode is described by {9, cj)) which 
correspond to the following fluctuations: (jf^ = {Sijd"^ — 
didj)(j) and 9 = {Sijd^ — didj)9ij. Its frequency is again 
Ldk = 0. The Hamiltonian for such a mode contains only 
9. So the quantum fluctuations satisfies 5(f> ^ 1 and 
S9 <C I /no. The second helicity mode is also gapped 
with a gap of order U2. The strong fluctuations S<j>''-^ ~ 
{Sijd^ — didj)4> ^ 1 lead to a gauge transformation 



(58) 



and the weak fluctuations 9 = {Sijd"^ — didj)9ij <^ I/jt-g 
lead to a constraint 







(59) 



18 



The key in our argument is the gapping of heUcity 
and ±1 modes. Clearly, in the L-type lattice model (60), 
those modes are strongly fluctuating and strongly inter- 
acting modes. Those modes also have a very narrow band 
width and arise from compact degree freedom on lattice. 
The strongly interacting modes with flat bands are in 
general gapped. The gapless helicity ±2 modes are very 
classical in the large uq limit. They should survive the 
gapping of the helicity and ±1 modes. 

VII. THE N-TYPE LATTICE MODEL 



• Giving non-trivial dynamics to states in the low en- 
ergy Hilbert space through Hj and Hg . 

In this section, we will try to construct a lattice model 
whose low energy effective theory is the N-type model 
(33) with the constraints cqn. (13) and eqn. (21).^^ Just 
like the L-type lattice model, the Hamiltonian of the N- 
type lattice model contains Hu (39) which imposes the 
constraints eqn. (13) and eqn. (21). Other terms of the 
Hamiltonian can be obtained by putting the continuum 
Hamiltonian (32) on lattice. 

A. Putting flijii*-' on lattice 

If wc introduce D\, = e^^^diOjk, we find that OijW^ ^ 
D\D'l. On lattice D[ can be chosen to have the following 
form 

Dy = dyUzy — dzCLyy ^ 2 ' 
= dyttzz - d^Oyz - , 

where Here Pj{i) are defined as 

P^(^) = Oz.ii + y + f + f ) + Ozx{i + f + f ) 

- e^yii + 2; + 1 + 1) - e,y{i + f + 1), 
;3^(i) = -0.,(i + | + |)-0,,(i + |-|) 

+ 29yy{i + Z)+2eyy{i), 

Pl{i) = -2e,,{i + y)-2e,,{i) 

+ 0yzii + ^ + ^) + Oyzii + ^-^). 

Other components are obtained by cycling xyz to yzx 
and zxy. So ^goijR^^ leads to the following contribution 
to the Hamiltonian 

i,a—x,y,z i,ab—xy,yz,zx 



In the compactificd model, the following term 

^^ = ^ E n-<^os[p:m 

i,a—x,y,z 
i.ab—xy,yz,zx 

reproduce H'g in the small p'^ limit. The operators in Hg 
are functions of Wl'' and Wg'^. So Hg acts within the 
Hilbert space of compactified model. 

B. Putting £'^£'^ - on lattice 

The term ^ J[£^'-'f — ^£"] can be put on lattice easily. 
It leads to the following contribution to the Hamiltonian 

H'=J y L-\i + ^ + ^)L-\i + ^ + ^) 
J ^ ^ 2 2^ ^ 2 2^ 

i ,ab—xy ,yz ,zx 
i.a—x^y^z 

In the compactificd model, the following term 
Hj^ncJ E {1 -cos[27rL°°(i)/nG]} 

%.a—x,y,z 

-f2nGjE E {I ~cos[2t: L''\i)/nG\} 

i ah—xy,yz.zx 

-inGjE{l-cos[2^ E 

i a—x,y,z 

reproduce Hj in the small L"^^ /nc limit. 



C. The N-type lattice model 

Putting everything together, we obtain the following 
lattice Hamiltonian 

H = Hu + Hj + Hg (60) 

which defines the N-type lattice spin model. We will 
assume U ^ J ^ g in the N-type lattice model. The 
corresponding phase-space Lagrangian is given by 

«G b.^^.. a b. 

L = ^ E ri^+^ + -^)^oO.>,i^ + ^ + -^) 

i,ab—xy,yz,zx 

+ if E r^{i)doeaaii)^Hi<l>^^,e,,) (61) 

i,a—x,y,z 

where c/)^^ = 2tt L^^ /uq and H{(j)^^,9ij) is given by (60) 
with L*^ replaced by nc'f)''-' /2tt. 
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D. Low energy collective modes — 
spin- wave / quantum- freeze approach 

• In a semiclassical approach to the N-type model, we 
can expand the classical action to the quadratic or- 
der and use the spin-wave/quantum-freeze approach 
to obtain the low energy dynamics. 

The motivation of constructing the quantum spin 
model (60) or (61) is to obtain a lattice model whose 
low lying collective modes are gravitational waves with 
a linear dispersion. After we obtain the model (60), we 
throw always the motivation and ask "do the low energy 
excitations of the spin model (60) really correspond to 
gravitons?" 

The spin model (60) is a complicated and strongly in- 
teracting quantum system. It seems impossible to ob- 
tain and to understand the dynamics of its low energy 
excitations. On the other hand, the rotor model (11) 
is also a complicated and strongly interacting quantum 
system. Many different approaches, such as equation- 
of-motion approach, coherent state approach, and spin- 
wave/quantum- freeze approach, are developed to under- 
stand the low energy dynamics of (Il).8.i8,i9,37 

In this section, we will use the spin- wave/quantum- 
freeze approach developed in section III D to understand 
the low energy excitations of the spin system described by 
the Hamiltonian (60) in large limit with g, J, Ui and 



U2 are of the same order. We would like to remark that 
the spin- wave/quantum- freeze approach has not been 
tested fully and may not be as reliable as the cquation- 
of-motion approach used in section A. We note that for 
the N-type lattice model, Hj and Hg do not act within 
the zero-energy subspace of Hjj (since [Hjj, Hj] ^ 0, and 
[i?(7,i/g] 7^ 0). The equation-of- motion approach cannot 
be used for the N-type lattice model. This is why we 
will use spin- wave/quantum- freeze approach to analyze 
the low energy excitations of the N-type lattice model. 
The results obtained in this section may not be reliable 
and need to be confirmed through other more reliable 
methods. 

In the spin- wave/quantum- freeze approach, we first as- 
sume the quantum fluctuations of Qij and to be much 
larger than their discreteness A6'ij ~ ^4>ii ~ ^/nc- (We 
will check this assumption later for self consistency.) In 
this case, we can treat 9ij and <f>ij as classical fields and 
use the classical phase-space Lagrangian (61) to study 
the low energy collective modes. 



1. Classical spin waves 

Let us also assume that the quantum fluctuations of 9ij 
and 4>ij are much smaller then 1 and expand the phase- 
space Lagrangian to quadratic order in Oij and (pij. In 
fc-space, the quadratic Lagrangian has a form; 
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where Cx.y.z and Cxy 



are defined as, for example Cx 



COS ^ COS ^ 



After solving the equation of motions obtained from 
the quadratic Lagrangian, we can find the dispersions of 
collective modes of our bosonic model (60). There are to- 
tal of six collective modes. We find four of them have zero 
frequency for all fe, and two modes have a linear disper- 
sion relation near k = (tTjTT, tt). Near k = (tTjTt, tt), the 
dynamics of the six modes are described by the following 
continuum field theory: 



C^na 



I 2^' 



J 



m\2 



ncUi 



M\2 



2 

nGU2 



ii\2 



{R") 



(62) 



where i?*-' = e™''e^'"9m9za„A;. The continuum fields 
abcix) are given by — (— l)*i0'"^(i + | + |) for & ^ c and 
by {—iyO^'^{i) for b = c, where b,c = x, y, z. The contin- 
uum fields (l)°-''{x) are given by + f + |) for 
a^b and by (-l)V°''(i) for a = b. 

2. The classical heUcity ±2 modes 

The two modes that have a linear dispersion relation 
are the two helicity ±2 modes. The dispersion relation 

is given by ojk ~ We note that the U2 and 

Ui terms decouple from the helicity ±2 modes, since the 
helicity ±2 modes satisfy the constraints dicf)^^ ~ R" =- 0. 
Thus the dynamics of the helicity ±2 modes does not 
depend on C/2 and U\. 

To see if the result cjfc \J^\k\ can be trusted, we 

need to analyze the quantum fluctuations of Oy- and (jf '-i . 
We find that for large uq, the quantum fluctuations of the 
helicity ±2 modes is of order Scj)^^ ,Saij ^ \J\JnQ (note 
that J and g are of the same order). So the fluctuations 
of (^^^ and ay satisfy l/ric ^ ^'i>^^ ^ 1, l/"-G ^ ^o^ij 
1 and the semiclassical approximation is valid for the 



helicity ±2 modes. In this case, the result Uk ~ ^^J gJ\k\ 
can be trusted. The spin model (60) has gapless gravitons 
as its low energy excitations. 



that those modes are quantum frozen and are gapped. 
Since discreteness of (j)^^ is of order l/nc, the gap of the 
helicity and ±1 modes is of order Ui/nc- The strong 
fluctuations Saij — diOj + djOi 3> 1 lead to gauge trans- 
formations 

+ d,e, + d,e, (63) 

and the weak fluctuations 0* <C 1/ nc lead to constraints 







(64) 



The second helicity mode is described by (6*, 0) which 
correspond to the following fluctuations: (f)^^ = {Sijd"^ — 
didj)(j) and 6 = {Sijd^ — didj)aij. Its frequency is again 
ujk = 0. The Hamiltonian for such a mode contains only 
6. So the quantum fluctuations satisfies 50 3> 1 and 
S9 <C l/nc- The second helicity mode is also gapped 
with a gap of order f/2. The strong fluctuations S<j>''-^ ~ 
{Sijd^ — didj)4> ^ 1 lead to a gauge transformation 

^ cj)'' + {S,jd^ - d.,d,)(P (65) 

and the weak fluctuations 9 = {Sijd'^ — didj)aij <C I/'t.g 
lead to a constraint 



(%c>2 _ d,d,)a,j = 



(66) 



The Lagrangian (62) equipped with the gauge trans- 
formations (63,65) and the constraints (64,66) is nothing 
but the linearized Einstein Lagrangian of gravity, where 
aij ~ gij — 5ij represents the fluctuations of the metric 
tenor gij around the flat space. So the linearized Einstein 
gravity may emerge from the quantum model (60) in the 
large uq limit. The qubit model (60) may be quantum 
theory of gravity at the linear order. 

The key in our argument is the gapping of helicity 
and ±1 modes. Clearly, in the N-type lattice model (60), 
those modes are strongly fluctuating and strongly inter- 
acting modes. Those modes also have a very narrow band 
width and arise from compact degree freedom on lattice. 
The strongly interacting modes with flat bands are in 
general gapped. The gapless helicity ±2 modes are very 
classical in the large uq limit. They should survive the 
gapping of the helicity and ±1 modes. 



3. The quantum freeze and the gapping of helicity and ±1 

modes 

The helicity ±1 modes and one of the helicity mode 
are described by (0^, (pi) which correspond to fluctuations 
of the following form: aij = diOj -\- djOi and 0* = djcfP'^ . 
Their frequency ujk = 0. For such modes, the Hamilto- 
nian only contains Thus the quantum fluctuations 
satisfies (50' <C l/?^G and 59 i 1. So the semiclassical 
approximation is not valid and the result Wfc = cannot 
be trusted. Using the similar argument used in the emer- 
gence of J7(l) gauge bosons in section HID, we conclude 



E. Difficulties with the N-type lattice model 

• Since [Hu^Hj] ^ 0, and [Hii^Hg] ^ 0, the low en- 
ergy dynamics of the N-type model is hard to obtain 
reliably. The spin-wave/quantum-freeze approach 
is an uncontrolled approximation, since the terms 
beyond the quadratic order are not unimportant. 

We have seen that the helicity ±2 modes are very clas- 
sical with small quantum fiuctuations while the helicity 
and ±1 modes have strong quantum fluctuations. This 
is why the helicity ±2 modes remain gapless while the 
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helicity and ±1 modes are gapped. Also, at quadratic 
level, the helicity ±2 modes and the helicity 0, ±1 modes 
decouple. The strong quantum fluctuations in the 0, ±1 
modes will not affect the dynamics of the ±2 modes at 
the quadratic level. 

However, due to the strong quantum fluctuations in the 
helicity 0, ±1 modes, the quadratic expansion for those 
modes is no longer valid. Beyond the quadratic level, the 
helicity ±2 modes and the helicity 0, ±1 modes are cou- 
pled in the N-typc lattice model. In this case, the strong 
quantum fluctuations in the 0, ±1 modes may change the 
dynamics of the ±2 modes beyond the quadratic level. 
This is why we cannot obtain the low energy dynamics 
of the N-type lattice model reliably and the a; |fe| dis- 
persion for the helicity ±2 modes in the N-type lattice 
model is not a reliable result. 

However, such problem does not appear in the L-type 
lattice model. The helicity ±2 modes and the helicity 
0, ±1 modes decouple even beyond the quadratic level, 
since Hu commutes with Hj + Hg in the L-typc lattice 
model (54). So the u ~ \k\^ dispersion for the helicity 
±2 modes in the L-type lattice model is a reliable result. 



VIII. TOPOLOGICAL DEFECTS 



In section VC, we have constructed the low energy 
Hilbert space of our lattice model in the large Ui^2 limit. 
We can also construct high energy states that do not 
satisfy the scaler or vector constraints (41) but remain to 
be eigenstates of S{i) and Q(i, i + a). These states can 
be naturally regarded as scalar defect and vector defect 
as discussed in the continuum limit in section IV. 



with 



For example: 



\n=0 



is a scalar defect with 



i27rn- 



SiiT |0) 



(67) 



Til 

S{j)\m,i) = Si j exp 27ri — 
m = 0,l,2---riG - 1 



TO, I) 



(68) 



The vector defect can be constructed as, for example: 



qx.lxy = 0, l,2---nG - 1 



• |0) 



Q{3J + a)\q^,i) 

= [5i^j5x.a. + Si-ajSx,a) exp 



2TTi 



no 



\qx,i) (70) 



QUJ + a)\lxy,i) 

.1 



exp 



27ri- 



\lxy, i) 



(71) 



Here the operator (^X]nSo*^^P 



(69) 



notice the vector defect can only be created or annihilated 
in pairs. Thus the vector defect in our lattice model has 
a fractionalizcd nature. 

' S{iY ) and 

(M^e^^(i))'=" , [Wg^ii + % + 1)'°=" can be viewed as scalar 
and vector defects creating operators. The quantum 
number m,qa,lab divide the Hilbert space into differ- 
ent sectors, our low energy Hilbert space is defined as 
m ^ Qa ^ lab = 0. 

The existence of (quantized) scalar and vector defects 
have very deep relationship with the existence of helicity 
±2 modes in low energy Hilbert space. In the continuum 
limit, it is easy to see these defects are all topological de- 
fects and could not be removed by local operations and 
protected by energy scale Ui and U2- Meanwhile, this 
energy scales "freeze" the low energy quantum fluctua- 
tions of helicity and ±1 modes. At this stage, we may 
conclude the helicity ±2 modes in our lattice model is 
topologically stable and the existence of scalar and vec- 
tor defects reflects the special quantum order in such a 
system. 



IX. CONCLUSION AND DISCUSSION 

In this paper, we studied two qubit models, the L-type 
model and the N-typc model. The low energy dynam- 
ics of the L-type model can be calculated reliably and 
we find that the helicity ±2 modes are the only gapless 
excitations with uj ^ |fc|^ dispersion. The low energy 
dynamics of the N-type lattice model cannot be calcu- 
lated reliably. However, the spin- wave/quantum freeze 
approach suggest that the helicity ±2 modes may be the 
only gapless excitations with w ^ |fc| dispersion in the 
N-type model. Those gapless helicity ±2 excitations can 
be viewed as emergent gravitons. 

Both qubit models are constructed by putting a sym- 
metric tensor field on lattice, and by discretizing and 
compactifying each component of the tensor field. To 
understand the low energy dynamics of the two model, 
we have used a semiclassical approach, where we treat 
the tensor field as continuous and uncompactified. We 
then expand the lattice Lagrangian of the two models to 
quadratic order. The resulting free theories can be an- 
alyzed easily. We find that both models contain gapless 
helicity ±2 modes. The analysis of the quantum fiuctua- 
tions for such modes indicates that the helicity ±2 modes 
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are very classical and the semiclassical approach is self 
consistent. We also find that both models contain gap- 
less helicity 0, ±1 modes. The analysis of the quantum 
fluctuations for those modes indicates that the helicity 
0, ±1 modes are extremely quantum and the semiclassi- 
cal approach is not self consistent. The non-linear effects 
of the helicity 0, ±1 modes need to be considered. 

However, the non-linear effects of the helicity 0, ±1 
modes have very different forms in the L-typc and the 
N-type model. In the L-type model, the strong fluctuat- 
ing helicity 0, ±1 modes decouple from the semiclassical 
helicity ±2 modes even at non-linear level. After tak- 
ing into account the discreteness and the compactness of 
the tensor field, we can show that non-linear effect gaps 
the helicity 0, ±1 modes. Such a gapping process does 
not affect the helicity ±2 modes. So helicity ±2 modes 
remain gapless with cj ~ dispersion. 

In the N-type model, the strong fiuctuating helicity 
0, ±1 modes start to couple to the semiclassical helicity 
±2 modes at non-linear level. After taking into account 
the discreteness and the compactness of the tensor field, 
we can argue that non- linear effect gaps the helicity 0, ±1 
modes. But such a gapping process may affect the low 
energy dynamics of the helicity ±2 modes. This is the 
reason why our result for the N-type model is not reliable. 

The situation for the emergence of the uj k grav- 
ity from the N-type model is very similar to the emer- 
gence of Chern-Simons theory from local qubit models in 
1-I-2D. Although we have no reliable analytic calculation 
to show that the Chern-Simons theory can emerge from 
qubit models, numerical calculations have convincingly 
established the emergence of Chern-Simons theory from 
certain qubit models. So similarly, it is possible that the 
N-type model can produce the linearized Einstein grav- 
ity, even though we can not show it reliably. Numerical 
calculations are needed to show it convincingly. 

We have seen that the gapping of the helicity mode in 
the rotor model (11) leads to an emergence of J7(l) gauge 
structure at low energies. The emergence of a gauge 
structure also represents a new kind of order - quantum 
order*''"' - in the ground state. In Ref. 41, it was shown 
that the emergent U{1) gauge invariance, and hence the 
quantum order, is robust against any local perturbations 
of the rotor model. Thus the gaplessness of the emergent 
photon is protected by the quantum order. 

Similarly, the gapping of the two helicity modes and 
the helicity ±1 modes in the bosonic model (or qubit 
model) (54) and/or (60) leads to an emergent gauge 
invariance of the linearized coordinate transformation. 
This indicates that the ground state of the qubit models 
contain a new kind of quantum order (which is similar 
but different from those associated with emergent ordi- 
nary gauge invariance of internal degrees of freedom). 
We expect such an emergent linearized diffeomorphism 
invariance to be robust against any local perturbations 
of the bosonic model. Thus the gaplessness of the emer- 
gent helicity ±2 modes are protected by the quantum 
order and are robust against any local perturbations. 



Let us discuss the robustness of the uj ^ k^ dispersion 
of the L-type model in more detail. As stressed before, 
uj ^ k^ dispersion is the direct result of the local gauge 
invariance under eqn. (15) and eqn. (23). To be more 
precise, we know that the L-type model is described by 
the following effective Lagrangian 

C = r^9oa,j - - ^R'^R'^... (72) 

Its Hamiltonian density is given by 7^ = jCpj + 
^_R*-'i?*-' + ..., where ... represent higher order and/or non- 
quadratic terms. Such an effective Hamiltonian density 
for the L-type model is directly invariant under the gauge 
transformations (15) and (23), which leads to the a; ^ fc^ 
dispersion. Now, we would like to ask, if we integrate 
out high energy fluctuations in the L-type model, can we 
generate to following form of Hamiltonian density 

J7^ = y[(r^)'-i(rf] + ^a.,i?'^- 

= + (73) 

that appears in the N-type model? Here D'j^ = e^^^diajf^. 
If such Hamiltonian density is generated, it will modify 
the low energy dispersion from uj k^ to cu ^ k. 

Here we would like to stress that although dH = 
J cPx 5% is invariant under the gauge transformations 
(15) and (23) (with the constraints (13) and (21)), the 
Hamiltonian density Sl-L is not invariant under the gauge 
transformations. 5% changes by a total derivative. Since 
% is locally gauge invariant, within a perturbative calcu- 
lation, it cannot generate 5% that is only gauge invariant 
up to a surface term. So the to ^ k^ dispersion of the L- 
typc model is robust against all perturbative corrections. 
At the moment, we do not know any non-perturbative 
effects that can generate the 5T-L term within the L-type 
model. 

We further find the existence of scalar and vector de- 
fects to have intrinsic relationship with the gapping of 
and ±1 modes. At this stage, we may interpret the scalar 
defect as mass and vector defects as momentum. They 
are all quantized topological defects. However, there's no 
direct evidence for mass quantization experimentally, we 
attribute this fact to two possible reasons: a) the mass 
quanta is too small to be detect experimentally; (of or- 
der U2/TLQ in the large tlq limit) b) the appearance of 
cj ~ |fc| gravitons may remove the mass quantization for 
low lying excitations. To solve these problems demands 
much deeper understanding of the origin of the linear 
dispersion relations and will be our future efforts. 

This research is supported by the Foundational Ques- 
tions Institute (FQXi) and NSF Grant DMR-0706078. 
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Appendix A: Low energy collective modes for L-type 
lattice model — equation-of-motion approach 

In this section, we will use a different approach - 
equation-of-motion approach to study the low energy 
modes of the L-type spin model (54) . The idea behind the 
equation-of-motion approach is to choose operator that 
act within the low energy constrained subspace, and find 
the equation of motion for the averages of those opera- 
tors. The derived equation of motion will describe the 
classical motion of the low energy collective modes, pro- 



vided that the operators are chosen properly. For our 
case we will choose the gauge invariant operator TZab and 
T§ to construct the equation of the motion. 

In the large uq limit, we can express exp(27ri/nG') as 
1 -f- 27ri /uq and the commutating relation becomes 

Wab^ab _ y^ab^ab ^ WS^Wt (Al) 

no 

notice and TZab commute with the Hu term, so the 
equation of motions can be simplified as, for example: 



dt{ny.{i + I + f )) = -m.T^y.ii + f + f )]) = -i{[H^,nyS + I + |)]) = 

(~T|(i + y) - r;(i -y)+ T^ii + z) + T^{i - z) - T^{i - x + y) - T^{i - x - y) + T^ii - x + z) 

+ T^{i ~x-z)- 5Ty{i) - bTy{i -x)- Ty{i + x)- Ty{i - 2x) - Ty{i + z)- Ty(i -z)- Ty{i -x + z) 

- Ty{i -x-z) + 5T^{i) + 5T^{i -x)+ T^_{i + x) + T^(i - 2x) + T^(i + y) + T--{i - y) + T^{i - x + y) 

+ T:{i-x-y) + ih4T^(i) - 4Ty{i + y) - + z) - Ty{i - z) - Ty{i + x) - Ty{i - x) - Ty(i + z + y) 

- Ty{i -z + y)-Ty{i + x + y)- Ty{i -x + y)- AT^ii) - AT^{i + y) - T^{i + z) - T^{i - z) - T^ii + x) 

- T^{i -x)- T^{i + z + y)- T^{i -z + y)-T^{i + x + y)- T^{i - x + y) + 4T^(i) + 4T^{i + z) + T^{i + y) 
+ T^ii - y) + (i + x)+ T^{i -x)+ T^{i + y + z) + T^{i - y + z) + {i + x + z) + T^{i - x + z) + 4:T^{i) 
+ AT^{i + z)+ T^{i + y)+ T~{i - y) + T^{i + x)+ T^{i - x) + T^,{i + y + z)+ T^{i - y + z) +T^{i + x + z) 

+ T^{i -x + z)]- h.c.){ny,{i + I + f )) (A2) 



dt{n..{i)) = -^{[H, n,,{i)]) = -^{[H^,^Z.,{^)]) = - — {-2[T^{i) + T^{i - x) + T^. {i - y) + r;(i - z) 
+ T|(i -x-z)+T^{i-y-z)+ T^ii - x - y) + T^{i - x - y - z)] + 2[Ty{i) + Ty{i - x) + Ty{i - y) 

+ Tyii - z) + Tyii - X - z) + Tyii - y - z) + Ty{i - X - y) + Tyii - X - y - z)] + ^[Tyii + y) + Tyii - y) 

- Ty\% + x)- Ty\i -x) + Ty(i -z + y)+ Ty{i -z-y)- Ty{i -z + x)- Ty{i -z-x)+ T^ii + y) 
+ Ty^ii-y)-T^{i + x)-Ty^i-x)+T^'{i-z + y)+T^{i-z-y)-T^{i-z + x)-T^{i-z-x) 
+ 5T^{i) + bT^ii -x)+ T^{i + x)+ T^a - 2x) + T^{i + y) + T^^(i - y) + T^{i + y - x) + T^{i - y - x) 
+ 5Ty{i) + 5Ty{i -x)+ Ty{i + x)+ Ty{i - 2x) + Ty{i + y) + Ty{i -y)+ Ty{i + y-x)+ Ty{i -y-x) 

- bT^ii) - bT^ii - y) - T^ii + y) - T^ii - 2y) - T^ii + x) - T^{i - x 

- bT^ii) - bT^Ji - y) - T^a + y) - T^{i - 2y) - T.^{i + x) - T^{i - x) 
-h.c.){nUi)) (A3) 



- Tf (i + x-y)- T^{i -x-y) 
T^{i + x-y)-T^{i-x-y)] 
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dt{T^{i)) = -z{[H,T^ii)]) = -^{[H^,T^{^)]) = -|-(T:(i))x 

{2[nyy{i) + TZyy{i + X)+ 7^,y,y(^ + y)+ Uyyil + Z) 

+ Uyyii + x + y)+ TZyy{i + y + z) + TZyy{i + z + x) + TZyy{i + x + y + z)] 

- 2[7^,,(^) + 7^,,(^ + x) + n,,{i + y)+ n,,{i + z) 

+ n,,{i + x + y)+ n,,{i + y + z)+Ji:,^{i + z + x)+Ti,,,{i + x + y + z)] 

+ bUxyii) + 5TZxy(i + z)+ TZ^.y{i + 2z) + TZ^.y{i - z) + Uxyii + y) + Hxyii - y) + 'R-xyii + z + y) 

+ Tlxyii + z-y)+ TZyzii + y) + 7^y^(^ - y) - TZyz{i + z) - 7^y^(^ - z) + 7^y^(^ + x + y) + 7^^,^(^ + x - y) 

- 'R-yzii + X + z)- TZyz{i + x- z)- bTZzxii) - 57^^J,(^ + y) - TZzx{i + 2y) - TZzx{i - y) - Tizxii + z) 

- n,x{i -z)-n,x{i + y + z)- n,x{i + y-z)- h.c.) (A4) 



d,{T^{i)) = r;(i)]) = -^([i/^^T;(^)]) = V(r,-(i))x 

{-bTlxx{i) - ^TZxxii + z) - TZxxii + 2z) - Tlxx{i - z) - Tlxx{i + y) - Tlxx{i - y) - 'Tixx{i z + y) - TZxxii + z - y) 

+ 57iyy{i) + 57^yy(i + z) + TZyy{i + 2z) + Tiyy{i - Z) + TZyy{i + X) + 7lyy{i - x) + TZyy{i + Z + x) + TZyy{i + Z - X) 

+ Tlzz{i + y) + Uzzii - y) - nzz{i + x)^ TZzzii - x) + TZzzii + z + y) 
+ Tlzz {i + z~y) - TZzz {i + z + x) - TZzz {i + z - x) 

+ 47^J,^(^) + ATZyzii - y) + 7^y^(^ + a;) + TZy^ii - x) + 7^y^(^ - y + x) + TZy^ii - y - x) 
+ Uy^Xi + z) + TZyzii - z)+ 7^y^(^ -y + z)+ TZy^ii - y - z) 

- ^T^zx{i) - 47^^^(^ -x)- n.,^,{i + y) - 7^^2;(^ - y) - Tl^.^{i - x + y) - TZ,_^{i - x - y) 

- T^zxii + z)- n^x{i - z)- n^^^ii -x + z)- n^r,{i -X- z) ~ h.c.) (A5) 



dt{T:ii)) = Tf(i)]) = -^([i^^^T;(^)]) = ^/(Tf (i))x 

{^TZxxii) + 5TZ^x{i + y) + TZxxii + 2y) + 7^J,2;(^ - y) + 7^J,:^;(^ + z) + 7^J,^(^ - z) + TZxxii + y + z) + TZxxii + y - z) 

+ Tlyy{i + X)+ TZyy{i ^X)- Tiyy{i + Z) - 'Ryy{i - Z) 

+ Uyy{i + y + x)+ Uyy{i + y ^ x) ^ TZyy{i + y + z)- TZyy{i + y - z) 

- 57^^^(^) - 57^^^(^ + y) - 7^^^(^ + 2y) - TZzzii - y) - 7ez^(i + x) - 7^^^(^ - a;) - 7^^^(^ + y + a;) - 7^^^(^ + y - a;) 
+ 47^:l,J,(^) + ATZxyii - x) + TZxy{i + y) + 7^3;y(^ - y) + "^ajyl* - a; + y) + 7^l,J,(^ - a; - y) 

+ TZxyii + z) + 7?.2:y(i - z) + TZxy{i - x + z) + TZxy{i - x - z) 

- ^TZyzii) - 4:TZyz{i ~ z) - TZyz{i + y) - TZyz{i - y) - TZyzii - z + y) - 7^y^(^ - z - y) 

- 7^yz(i + a;) - TZyz^i - x) ~ TZyz{i - z + x) - TZyz{i - z - x) - h.c.) (A6) 

I 



other components of TZab and can be obtained by sim- 
ply cycling xyz to yzx and zxy. 

In the large nc limit, the fluctuations A0a6 ^ 
l/nc is much less than 1. So {T^) ~ (Teab) ~ 1 (notice 
and T^-ab are unitary operators). One may express 



as: {TZab) = e'®-'' , where Qab ^ is a small fluctuations 
around 0. This allow us use the approximation sin Qab ^ 
Qab and linearized the above equation of motions in oj, fc 
space. Introducing 
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where c^; = 2 cos ^ , Cy = 2 cos ^ , = 2 cos ^ 
write the hnearized equation of motions as 



Qyy (^) 

V e,,(fe) / 



e,,(fe) 

0a;i^ (fc) 



By solving the hnearized equation, we find that four of 
the six modes have zero frequency. Those four modes 
have strong quantum fluctuations and wiU be gapped 



once we take into account the compactness and discrete- 
ness of the fields Qab- The only two low lying dynamical 
modes are helicity ± 2 modes with the dispersion rela- 
tions: 



2Tr\ 



h 



(A7) 



This allows us to conclude that the low energy effective 
theory of the lattice bosonic model (54) is eqn. (28) with 
the constraints eqn. (13) and eqn. (21). The low energy 
excitations are helicity ±2 modes with a; cx fc"^ dispersion. 
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A lattice bosonic model and a lattice qubit model (ie a 

lattice spin model) are equivalent. So in this paper the 

two terms, bosons and qubits, are equivalent and are used 

interchangeably. 

In this paper, we define gravitons as linearly dispersing 
gapless helicity ±2 excitations. 

Here, quantum theory is defined as a theory whose states 
are state vectors in a Hilbert space, and observables and 
Hamiltonian are Hermitian linear operators. So in quan- 
tum theory, by definition, there exists an absolute time. 
[We note that path integral formulation can some times go 
beyond quantum mechanics (say when time is discretized).] 



